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ABSTRACT 

We  have  inereased  the  spectral  resolution  and  accuracy  of  the  absolute  infrared  spectral  flux  for  a  subset  of  stellar 
standards  created  by  Cohen  and  his  collaborators  in  previous  papers  in  this  series.  We  combined  the  moderate- 
resolution  (A/AA  ~  400)  spectroscopy  obtained  on  infrared  standard  stars  by  the  Short  Wavelength  Spectrometer  on 
the  Infrared  Space  Observatory  with  high-quality  photometry  that  is  tied  to  our  recent  absolute  calibration  from  the 
Midcourse  Space  Experiment.  New  spectra  were  created  for  33  stars,  9  of  which  are  Cohen  et  al.  secondary  standards 
and  another  20  of  which  are  tertiary  standards  for  which  Cohen  adopted  spectral  templates.  The  effective  temper¬ 
atures  and  angular  diameters  for  the  stars  are  derived  from  the  absolute  spectra  and  compare  favorably  with  inde¬ 
pendent  measures  of  these  quantities  in  the  literature.  The  present  spectra  display  systematic  differences  with  those  of 
Cohen  and  colleagues  in  that  they  have  4%-7%  lower  fluxes  in  the  1-4  /rm  spectral  region.  Our  spectra  remove  the 
A-K  star  calibration  bias  recently  noted  in  the  calibration  of  the  Spitzer  Space  Telescope  camera. 

Key  words:  infrared:  stars  —  methods:  analytical  —  techniques:  spectroscopic 
Online  material:  extended  figure  set,  appendix 


1.  INTRODUCTION 

High-quality  space-based  spectroscopy  and  radiometry  were 
combined  to  improve  the  resolution  and  absolute  flux  accuracy 
of  the  infrared  spectra  for  a  subset  of  the  standard  calibration  stars 
created  by  the  pioneering  efforts  of  Cohen  and  his  collaborators. 
In  their  initial  paper,  Cohen  et  al.  (1992a,  Paper  I)^  defined  their 
primary  standards  as  Vega  (a  Lyr)  and  Sirius  (a  CMa).  They  es¬ 
tablished  the  zero-magnitude  flux  scale  by  extrapolating  the  ab¬ 
solute  flux  of  Vega  at  0.5556  pm,  recommended  by  Hayes  (1985), 
into  the  mid-infrared  using  the  calculated  spectral  flux  of  a  Kurucz 
(Teir  =  9400  K)  model  atmosphere.  Next,  they  adopted  a  Kurucz 
(Tesf  =  9850  K)  model  for  the  infrared  spectral  energy  distribu¬ 
tion  of  Sirius,  which  was  placed  on  an  absolute  scale  by  direct  in¬ 
frared  photometric  reference  to  Vega.  In  a  subsequent  series  of 
papers,  Cohen  et  al.  (1992b,  Paper  11;  1995,  Paper  IV;  1996b, 
Paper  VII)  created  absolutely  calibrated  composite  spectra  for  10 
bright  secondary  infrared  standards  by  seaming  together  mea¬ 
sured  spectral  fragments  across  the  near-  and  mid-infrared,  with 
each  fragment  referenced  to  the  primaiy  standards.  Cohen  et  al. 
(1996a,  Paper  VI)  adopted  a  model  atmosphere  for  the  spectmm 
of  q'  Cen  and  used  it  as  a  reference  to  create  secondary  composite 
spectra  for  two  additional  secondary  standards  in  the  southern 
hemisphere.  Ultimately,  Cohen  et  al.  (1999,  Paper  X)  populated 
the  sky  with  a  network  of  tertiaiy  standards  by  representing  their 
infrared  energy  distributions  with  “spectral  templates”  and  scal¬ 
ing  the  result  by  photometry  with  respect  to  either  a  primary  or 
secondary  standard.  The  template  for  a  given  star  was  taken  as  the 
composite  spectmm  of  the  secondary  standard  with  (nearly)  the 
same  spectral  type. 


‘  Institute  for  Scientific  Research,  Boston  College,  140  Commonwealth  Ave¬ 
nue,  Chestnut  Hill,  MA  02467;  charles.engelke.ctr@hanscom.af.mil. 

^  Air  Force  Research  Laboratory,  VSB,  29  Randolph  Road,  Hanscom  AFB, 

MA  01731;  kathleen.kraemer@hanscom,afmil,  steve.price@hanscom.afmil. 

^  We  refer  generally  to  the  stellar  calibration  and  standards  developed  by 
M.  Cohen  and  his  collaborators  as  “Coben  et  al.”  Specific  references  have  the 
year  of  publication  and  series  designation  for  the  first  instance  (e.g.,  1992a,  Paperl) 
and  the  series  designation  thereafter  (e.g..  Paper  I). 
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Price  et  al.  (2004)  assessed  the  accuracy  of  the  Cohen  et  al. 
calibration  using  the  results  of  the  extensive  series  of  calibration 
experiments  conducted  with  the  Midcourse  Space  Experiment 
{MSX).  MSX  obtained  a  direct  absolute  calibration  against  five 
emissive  reference  spheres  that  were  ejected  from  the  spacecraft 
at  different  times  during  the  mission  and  precisely  measured  the 
relative  photometry  between  the  Cohen  et  al.  primary  and  sec¬ 
ondary  standards.  To  summarize  their  results.  Price  et  al.  found 
the  following: 

1 .  The  MSX  experiments  confirmed  the  absolute  infrared  zero- 
magnitude  flux  scale  proposed  in  Paper  1  over  the  8-2 1  pm  range 
of  MSX  bands  A,  C,  D,  and  E.  The  absolute  calibration  against 
the  reference  spheres  agreed  with  the  Cohen  et  al.  absolute  zero- 
magnitude  infrared  fluxes  to  within  the  MSX  measurement  errors, 
which  were  less  than  the  1 .45%  uncertainty  Cohen  et  al.  ascribe  to 
their  zero-magnitude  fluxes. 

2.  The  MSX  relative  photometry  between  the  standard  stars 
agrees  well  on  average  with  that  derived  from  the  Cohen  et  al. 
composite  spectra  if  the  infrared  flux  from  Sirius  is  ~1  %  brighter 
than  the  flux  Cohen  et  al.  adopted  in  Paper  I.  The  formal  errors  and 
biases  of  the  MSX  photometry  for  the  standard  stars,  however,  are 
much  smaller  than  those  cited  by  Cohen  et  al.  (e.g.,  their  un¬ 
certainties  are  as  great  as  1 0%  in  the  center  of  the  4.3  pm  CO2 
atmospheric  band).  The  AfSYphotometry  can  therefore  be  used  to 
improve  the  accuracy  of  the  composites. 

The  low  (~0.1  ^m)  spectral  resolution  of  the  Cohen  et  al. 
composite  spectra  smooths  over  the  spectral  details  in  the  molec¬ 
ular  absorptions  in  the  atmospheres  of  the  standard  stars,  which 
can  introduce  errors  when  calibrating  narrow  spectral  bands  or 
spectral  measurements.  Sensors  that  use  very  narrow  spectral 
bands  require  photometric  standards  that  have  a  sufficiently  high 
spectral  resolution  to  accurately  reflect  the  spectral  features  that 
are  present.  The  higher  formal  accuracy  of  the  MSX  photometry 
can  be  used  to  reduce  these  uncertainties,  particularly  when  com¬ 
bined  with  higher  resolution  spectra  from  the  European  Space 
Agency’s  Infrared  Space  Observatory  (ISO;  Kessler  et  al.  2003) 
mission.  Here  we  describe  the  procedures  used  to  combine  the 
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TABLE  1 

Best  Secondary  Standards 


Tcs  6  Spectra" 

Star  Spectral  Type  (K)  (mas)  [to  X  (pm)]  Photometry*’ 

(1)  (2)  (3)  (4)  (5)  (6) 


PGem .  KOIII  4850  ±  5%  8.03  ±  3%  SEW,  SWS  9  Cen[9],  CVF  <5  Dra[16],  A[35]  MSX,  DIK,  H,  S,  T,  IRAS 

a  Boo .  K1.5  III  4350  ±  1.5%  21.06  ±  1%  SEW[2.0],  HW[2.4],  S[22],  A[35]  MSX,  DIR,  H,  S 

/?UMi .  K4III  4150  ±4%  10.00  ±  2%  S[17],  A[35]  DIR, /R4S 

a  Tau .  1C5  III  4050  ±  2%  20.75  ±  1%  SEW[2.0],  HW[2.4],  S[35]  MSX,  DIR,  H,  S 

7Dra .  K5  III  4030  ±  2%  10.17  ±  1%  HW[2.4],  S[26],  A[35]  MSA-,  DIR,  H,  S 

p  UMa .  MO  III  3900  ±  4%  8.45  ±  2%  HW[2.4],  S[17],  A[35]  DIR,  S,  T 

/3And .  MO  III  3900  ±  4%  13.65  ±  2%  SEW,  S[27],  A[35]  Am",  DIR,  H,  S 

QCet .  M2  III  3750  ±4%  12.94  ±  2%  SEW,  S[27],  A[35]  DIR, /7MS,  S 

7  Cm .  M4  III  3626  ±  2%  26.37  ±  1%  S[35]  MSX,  DIR 


Note. — Spectral  types  in  this  and  the  following  tables  are  from  Heras  et  al.  (2002)  and  references  therein. 

“  Spectral  segment.  (SEW)  Strecker  et  al.  (1979);  (HW)  Hinkle  et  al.  (1995),  Wallace  &  Hinkle  (1996,  1997);  (S)  SWS  (Sloan  et  al.  2003);  (CVF)  ISOCAM 
CVF  ( Engelke  et  al.  2004);  (A)  autoshape  (see  text).  The  number  in  brackets  is  the  wavelength  (in  microns)  at  which  that  segment  ends.  The  first  segment  starts  at 
1  pm  and  is  produced  with  autoshape.  Where  available,  the  Strecker  et  al.  (1979)  data  cover  1.2—2.36  pm;  otherwise,  autoshape  is  used.  The  SWS  spectra  all  start  at 
2.36  pm.  Adjacent  segments  start  where  the  previous  segment  ends. 

*’  Photometry  sources.  (MSX)  Price  et  al.  (2004);  (DIR)  Smith  et  al.  (2004);  (H)  Hammersley  et  al.  (1998);  (S)  Selby  et  al.  (1988);  {IRAS)  Beichman  et  al. 
(1988);  (T)  Tokunaga  (1984). 


MSX  photometry  with  ISO  spectra  to  produce  an  improved  set  of 
infrared  calibration  spectra.  Section  2  gives  an  overview  of  our 
results;  §  3  describes  the  data  processing;  §  4  presents  a  single¬ 
parameter  function  autoshape,  which  was  developed  to  fit  the 
measured  spectra  and  to  create  spectral  templates  for  stars  with 
good  photometry  but  no  spectra.  Section  5  compares  the  stellar 
parameters  derived  with  autoshape  to  those  in  the  literature;  §  6 
discusses  the  uncertainties  in  the  derivations;  and  §  7  summarizes 
our  results. 

2.  NEW  CALIBRATED  REFERENCE  SPECTRA 

The  model  spectral  energy  distributions  used  by  Cohen  et  al. 
for  Vega  and  Sirius  (Sirius  is  scaled  here  by  1 .01 ;  see  Price  et  al. 
2004  for  details)  are  adopted  for  this  analysis  in  order  to  preserve 
the  Cohen  et  al.  absolute  zero-magnitude  flux  scale  in  our  work, 
even  though  Short  Wavelength  Spectrometer  (SWS)  measure¬ 
ments  were  available  for  these  two  primary  standards.  The  SWS 
spectra  for  Vega  and  Sirius  are  noisy  and  plagued  by  detector 
nonlinearities  at  wavelengths  longer  than  about  7  )im.  Price  et  al. 
(2004)  obtained  an  absolute  direct  calibration  of  Sirius,  and  the 
photometry  in  the  present  analysis  is  referenced  to  this  standard. 
Although  they  also  directly  calibrated  the  infrared  flux  of  Vega,  it 
should  only  be  used  as  a  standard  at  A  :S  12  pm,  as  MSX  detected 
flux  excesses  from  the  known  debris  disk  at  longer  wavelengths. 
(In  Paper  II  Cohen  et  al.  recommended  limiting  the  use  of  a  Lyr 
below  15  pm  based  on  Infrared  Astronomical  Satellite  [/R/45] 
measurements.) 

We  have  combined  the  moderate-resolution  2.36-35  /rm  spec- 
tra(A/AA  ~  400)  obtained  by  the  SWS  ( Leech  etal.  2003)  on /SO 
with  the  best  photometry  available  on  Cohen  et  al.’s  secondary 
and  tertiary  standards,  which  we  reference  to  the  rescaled  absolute 
infrared  spectrum  of  Sirius,  to  create  absolutely  calibrated  moderate- 
resolution  spectra  for  33  stars,  29  of  which  are  in  the  Cohen  et  al. 
calibration  network.  The  highest  quality  photometry  was  from 
MSX,  the  Difiuse  Infrared  Background  Experiment  (DIRBE)  on 
the  Cosmic  Background  Explorer  (COBE),  and  Hammersley  et  al. 
(1998).  MSX  and  DIRBE  were  space-based  experiments,  and  there¬ 
fore,  atmospheric  effects  do  not  affect  their  measurements.  If 
needed,  measurements  from  IRAS,  as  well  as  the  photometric 
resources  listed  in  Appendix  A  of  Walker  &  Cohen  (1998), 
were  also  used. 


The  nine  secondary  standards  with  the  most  accurate  absolute 
infrared  fluxes  are  listed  in  Table  1,  while  24  secondary  stand¬ 
ards  with  greater  uncertainties  are  listed  Table  2.  The  stars  in 
Table  1  are  original  Cohen  et  al.  secondary  standards  from  which 
they  spawned  the  spectral  templates  for  the  stars  in  the  calibra¬ 
tion  network.  Eighteen  of  the  24  stars  in  Table  2  had  previously 
been  assigned  spectral  templates  in  Paper  X;  (3  Peg,  a  Cohen  et  al. 
secondary  standard  with  a  composite  spectrum,  is  relegated  to 
Table  2  because  it  is  a  small-amplitude  variable  star  in  the  infrared 
(Price  et  al.  2004).  The  derived  parameters  of  effective  tempera¬ 
ture  Teff  and  angular  diameter  9  and  their  associated  errors  are 
listed  in  the  tables  and  are  discussed  in  §  5.  References  for  the 
photometry  and  spectroscopy  used  to  create  the  improved  spec¬ 
tra  are  also  given  in  the  tables. 

The  absolute  spectral  energy  distributions  for  the  stars  in 
Table  1  were  constructed  using  precise  photometry  from  MSX 
and  DIRBE.  Because  the  stars  are  bright,  they  also  have  high 
signal-to-noise  ratio  (S/N)  SWS  spectra  over  much  of  the  wave¬ 
length  range.  Walker  &  Cohen  (2002)  created  a  template  spectrum 
for  one  of  these  stars,  /?  UMi,  that  they  then  used  as  the  rep¬ 
resentative  K4  III  template.  We  have  replaced  this  template  with 
an  SWS  spectrum  to  17  pm,  so  an  absolutely  calibrated  measured 
spectrum  is  now  available  for  this  star.  MSX  obtained  precision 
photometiy  on  p  Gem,  and  it  is  included  in  Table  1  even  though  it 
does  not  have  an  SWS  spectmm.  The  SWS  spectra  of  6  Cen  and 
6  Eri,  two  stars  of  similar  spectral  classification  to  P  Gem,  were 
averaged  and  adopted  as  the  composite  2.36-9  pm  spectram,  to 
which  was  seamed  the  9-16  pm  ISOCAM  circular  variable  filter 
(CVF)  spectrum  of  6  Dra  to  extend  the  measurements  to  16  pm. 
The  result  is  in  excellent  relative  agreement  with  the  high-quality 
MSX  photometry. 

Table  2  contains  secondary  standard  stars  that,  for  a  variety  of 
reasons,  have  larger  uncertainties  in  their  derived  absolute  fluxes 
compared  to  those  in  Table  1 .  The  increased  uncertainty  may  arise, 
for  example,  from  greater  uncertainty  in  the  photometiy.  For  in¬ 
stance,  Walker  et  al.  (2004)  estimate  a  6%-9%  uncertainty  in  the 
IRAS  12  and  25  pm  photometry,  which  is  about  10  times  larger 
than  that  for  the  MSX  photometry.  It  could  also  be  due  to  lower 
S/N  in  the  SWS  spectra,  which  limits  the  range  over  which  the 
spectra  can  be  used,  as  is  indicated  in  the  tables.  Two  of  the  orig¬ 
inal  Cohen  et  al.  secondary  standards  have  been  relegated  to  Table  2: 
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TABLE  2 


Additional  Secondary  Standards 


Star 

(1) 

Spectral  Type 
(2) 

/eff 

(K) 

(3) 

8“ 

(mas) 

(4) 

Spectra'’ 

[to  X  (pm)] 

(5) 

Photometry' 

(6) 

Cen . 

G2  V 

5870  ±  4% 

8.51  ±  2% 

P[2.36],  S[9],  A[35] 

B,  E,  Th 

/?  Dra . 

G2  II 

5100  ±  8% 

3.35  ±  5% 

SEW,  S[9],  A[35] 

DIR,  2M,  IRAS 

h  Dra . 

G9  III 

4950  ±  8% 

3.35  ±  5% 

S[10],  CVF[16],  A[35] 

HJRAS 

6  Eri . 

KO  IV 

4900 

2.48 

S[9],  A[35] 

B,  C,  VdB,  IRAS 

B  Cen . 

KO  III 

4800  ±  4% 

5.46  ±  2% 

S[9],  A[35] 

DIR,  B,  C,  VdB,  IRAS 

a  UMa . 

KO  Ilia 

4790 

6.68 

SEW,  S[9],  A[35] 

DIR,  Ken,  IRAS 

C  Dra . 

K2  III 

4570 

3.09 

S[10],  A[35] 

H,  S,  IRAS 

Q  Ari . 

K2  III 

4500  ±  4% 

6.89  ±  2% 

S[9],  A[35] 

MSX,  DIR,  A,  CX 

7  And . 

K3  Ilb 

4200  ±  4% 

7.96  ±  2% 

SEW,  S[10],  A[35] 

DIR,  A,  IRAS 

a  Tuc . 

K3  III 

4200  ±  4% 

6.19  ±  2% 

S[15),  A[35] 

DIR,  C,  IRAS 

X  Gru . 

K3  III 

4200 

2.82 

S[ll],  A[35] 

C.IRAS 

<T  Oph . 

K2  II 

4100  ±  8% 

3.52  ±  5% 

S[10],  A[35] 

DIR,  H,  IRAS 

6  Psc . 

K5  III 

4050 

3.75 

S[12],  A[35] 

DIR,  S,  VdB,  IRAS 

7  Phe . 

K4/5  III 

3950 

6.76 

S[10],  A[35] 

DIR,  IRAS 

H  Sco . 

K5  III 

3850  ±  8% 

4.90  ±  5% 

S[ll],  A[35] 

DIR,  C,  IRAS 

6  Oph . 

MI  III 

3850 

10.23 

S[13],  A[35] 

DIR,  C,  IRAS 

AE  Get . 

Ml  III 

3850 

5.19 

S[ll],  A[35] 

DIR,  F,  IRAS 

6  Vir . 

M3  III 

3660 

10.71 

S[ll],  A[35] 

DIR,  Ken,  IRAS 

p  Per . 

M4  n 

3540 

15.50 

S[19],  A[35] 

DIR,  A,  IRAS 

It  Aur . 

M3  II 

3500  ±  8% 

9.05  ±  5% 

S[10],  A[35] 

DIR,  K,  IRAS 

/3  Peg . 

M2.5  III 

3490  ±  8% 

17.88  ±  5% 

SEW,  S[35] 

MSX,  DIR,  H,  S 

/3  Gru . 

M5  III 

3480 

27.80 

S[26],  A[35] 

DIR,  Th,  IRAS 

GZ  Peg . 

M4  III  +  A2  V 

3450 

7.82 

S[20],  A[35] 

DIR,  Ker,  IRAS 

6^  Lyr . 

M4  II 

3300  ±  8% 

11.50  ±  5% 

S[27},  A[35] 

DIR,  IRAS 

‘  Except  where  noted,  the  uncertainty  in  Tij  is  ±6%,  and  that  in  6  is  ±3%. 

Spectral  segment.  (P)  Pickles  (1998);  (SEW)  Strecker  et  al.  (1979);  (HW)  Hinkle  et  al.  (1995),  Wallace  &  Hinkle  (1996,  1997);  (S)  SWS 
(Sloan  et  al,  2003);  (CVF)  ISOCAM  CVF  (Engelke  et  al.  2004);  (A)  autoshape  (see  text).  The  number  in  brackets  is  the  wavelength  (in  microns) 
at  which  that  segment  ends.  The  first  segment  starts  at  1  fim  and  is  produced  with  autoshape.  Where  available,  the  Strecker  et  al.  (1979)  data  cover 
1.2-2.36  /xm;  otherwise,  autoshape  is  used.  The  SWS  spectra  all  start  at  2.36  /rm.  Adjacent  segments  start  where  the  previous  segment  ends. 

'  Photometry  sources.  (A)  Alonso  et  al.  (2000);  (MSX)  Price  et  al.  (2004);  (DIR)  Smith  et  al.  (2004);  (F)  Feast  et  al.  (1990);  (Ken)  Kenyon 
(1988);  (Ker)  Kerschbaum  &  Hron  (1994);  (H)  Hammersley  et  al.  (1998);  (S)  Selby  et  al.  (1988);  {IRAS)  Beichman  et  al.  (1988);  (VdB)  Van  der 
Bliek  et  al.  (1996);  (Th)  Thomas  et  al.  (1973);  (B)  Bouchet  et  al.  (1989,  1991);  (C)  Carter  (1990,  1993);  (CX)  Paper  X;  (2M)  2MASS  (Stmtskie 
et  al.  2006), 


P  Peg  because  it  is  a  low-amplitude  variable  star  and  a'  Cen  be¬ 
cause  we  replaced  the  model  spectrum  from  Paper  VI  with  a  cal¬ 
ibrated  SWS  spectrum  that  has  good-quality  data  to  9  p;m.  The 
table  also  includes  <5  Dra.  Although  this  star  was  not  included  in 
the  Cohen  et  al.  calibration  network,  Cohen  created  a  calibrated 
template  spectrum  for  it  for  use  as  an  ISO  calibration  standard.'* 
The  9-16  fim  ISOCAM  CVF  (Blommaert  et  al.  2003;  Engelke 
et  al.  2004)  calibrated  spectmm  of  S  Dra  is  spliced  onto  the  shorter 
wavelength  SWS  spectral  segments.  The  CVF  spectmm  for  this 
star  has  a  much  higher  S/N  than  the  SWS  data  but  at  much  lower 
spectral  resolution.  The  6  Dra  spectmm  is  an  important  tie  point 
for  the  warm-temperature  end  in  the  analysis  of  the  trends  in  the 
stellar  continua  and  the  molecular  absorption  profiles.  To  anchor 
the  definition  of  these  trends  at  the  cool  end,  a  calibrated  spectmm 
for  P  Gm,  a  bright  southern  M5  III  star,  is  derived  even  though  this 
star  is  a  low-amplitude  variable,  p  Dra  was  not  in  the  Cohen  et  al. 
network  but  is  included  here  because  it  has  both  an  SWS  spectmm 
and  a  high-quality  near-infrared  spectmm  (Strecker  et  al.  1 979), 
and  it  provides  information  on  the  sparsely  populated  early  G 
spectral  classes  for  this  analysis.  Finally,  Walker  et  al.  (2004)  re¬ 
cently  noted  that  6^  Lyr  has  little  infrared  variability  and  should 
provide  reasonably  “accurate  and  reliable”  calibration;  because 
it  has  a  good-quality  SWS  spectmm,  it  too  is  included  in  Table  2. 


■*  See  http://www.iso.vilspa.esa.es/users/expl_lib/ISO/wwwcal/. 


3.  DATA  PROCESSING 
3.1.  Rationalizing  the  SWS  Spectrometry 

The /SO  SWS  obtained  over  1250  full-range  2.36-45  pm  spec¬ 
tra  on  more  than  900  objects.  About  50  of  these  spectra  are  of 
objects  most  suitable  as  calibration  references:  infrared  nonvar¬ 
iable  stars  without  excess  emission  from  circumstellar  dust.  The 
SWS  data  are  available  from  the  ISO  data  arehive  in  the  form  of 
the  288  individual  spectral  segments  that  were  obtained  during 
a  single  observation  (288  =  12  detectors  x  12  spectral  bands  x 
2  scan  directions).  The  archive  also  has  both  the  single  spectmm 
created  by  Sloan  et  al.  (2003)  from  the  288  spectral  fragments 
for  each  observation  and  the  separate  segments  they  created  for 
each  of  the  12  bands.  Sloan  et  al.  flat-fielded  the  measurements, 
removed  glitches  that  typically  arise  from  cosmic-ray  hits,  and 
then  smoothed  and  resampled  the  spectmm  onto  a  standard  uni¬ 
form  wavelength  grid. 

The  SWS  was  originally  calibrated  against  several  of  the  same 
Cohen  stellar  standards  with  composite  spectra,  such  as  7  Dra, 
that  are  upgraded  in  the  present  analysis.  A  mean  relative  spec¬ 
tral  response  for  the  spectrometer  was  determined  from  the  SWS 
observations  of  a  small  number  of  Cohen  et  al.  primaiy  and  sec¬ 
ondary  standards,  as  well  as  other  stars  with  model  atmosphere 
spectra.  The  relative  response  thus  derived  was  scaled  to  abso¬ 
lute  values  with  calibrated  photometry  obtained  by  the  ISO  pro¬ 
gram  for  that  purpose  (Hammersley  et  al.  1998;  Hammersley  & 
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Fig.  1. — Comparison  of  the  normalized  mid-infrared  spectra  of  a  Tau. 
(a)  Composite  spectmm  from  Paper  II.  (b)  SWS  spectrum  from  Sloan  et  al.  (2003) 
with  the  calibration  from  the  SWS  pipeline,  ver.  10. 1 .  (c)  Recalibrated  SWS  spec¬ 
trum  from  this  work.  All  three  spectra  have  been  normalized  with  an  Engelke 
function  at  TcH  =  4050  K  and  5  =  21.01  mas. 


5  10  15  20  25 

Wovetength  (^im) 


Fig.  2. — Comparison  of  the  observed  photometry  for  a  Tau  to  that  predicted 
from  the  spectra  shown  in  Fig.  1 ,  normalized  to  Sirius,  (a)  Percent  difference  for 
the  Cohen  et  al.  spectmm  shown  in  Fig.  la.  (b)  Percent  difference  for  the  SWS 
spectrum  shown  in  Fig.  16.  (c)  Percent  difference  for  the  final  recalibrated  SWS 
spectmm  shown  in  Fig.  Ic.  The  crosses  show  the  MSX photometry,  the  squares 
are  from  DIRBE,  and  the  triangles  are  from  Hammersley  et  al.  (1998). 


Jourdain  de  Muizon  2003).  The  SWS  team  estimated  that  the  rel¬ 
ative  spectral  response  has  an  overall  uncertainty  of  ~10%  but 
has  larger  uncertainties  of  ~20%-30%  at  the  longest  wavelengths 
(Leech  et  al.  2003). 

Sloan  et  al.  (2003)  spliced  the  12  SWS  spectral  segments  mea¬ 
sured  on  an  object  to  form  a  single  spectrum  by  assuming  that  the 
differences  between  segments  in  the  regions  of  overlap  were  either 
additive  or  multiplicative.  However,  noise  or  the  poor  responsiv- 
ity  in  the  overlap  regions,  as  well  as  incorrect  assumptions  about 
the  scaling,  could  cause  the  final  spectmm  to  depart  from  the  cor¬ 
rect  values,  occasionally  by  as  much  as  1 0%-20%  at  the  band 
edges.  Therefore,  for  a  given  star,  we  start  with  the  12  discrete 
spectral  segments  from  Sloan  et  al.  (2003).  These  segments  are  the 
average  of  all  the  “up”  and  “down”  spectra  fi-om  each  of  the  12 
detectors  in  the  sub-band,  which  improved  the  S/N  by  a  factor  of 
about  5  over  that  from  a  single  scan.  The  Sloan  et  al.  spectral  seg¬ 
ments  still  contain  residual  instrumental  artifacts,  due  to  factors 
such  as  residuals  left  by  the  nonlinear  response  correction  or  im¬ 
precise  background  subtraction,  which  are  identified  and  removed 
in  the  present  analysis.  Also,  the  spectra  become  increasingly  noisy 
beyond  8  /zm,  especially  for  the  fainter  stars  in  Table  2. 

The  stellar  flux  decreases  by  more  than  4  orders  of  magnitude 
across  the  2.36-45  fim  spectral  range  spanned  by  the  SWS.  The 
first  step  is  therefore  to  “flatten”  the  SWS  spectmm  by  dividing 
out  an  approximate  continuum.  We  approximate  the  infrared  stel¬ 
lar  continuum,  as  many  others  have  done  (e.g.,  Heras  et  al.  2002; 
Onaka  et  al.  2003;  Marengo  et  al.  1999),  with  the  Engelke  func¬ 
tion  (Engelke  1992).  This  function  is  based  on  the  semiempirical 


solar  spectmm  of  Vemazza  et  al.  (1976)  and  is  scaled  to  the  ef¬ 
fective  temperature  and  angular  diameter  9  for  a  given  star: 

ff.(>i.nff)  =  11,91  onr^ 

j  [ _ 19,500 _ 

_A7'eff(l-f  79450/27’eff)“  '®^ 

(1) 

where  H  =  n{9/2f  in  steradians.  This  function  is  a  good  repre¬ 
sentation  of  the  stellar  continuum  at  wavelengths  A  >  2  fim  for 
stars  of  solar  temperature  or  lower.  Price  et  al.  (2004)  found  that 
the  Engelke  continua  reproduced  the  MSX  photometry  on  the 
calibration  stars  to  a  few  percent. 

This  normalization  conveniently  scales  the  entire  spectmm 
to  a  single  well-resolved  plot  and  increases  the  visibility  of  the 
small-scale  local  trends  and  other  discrepancies  that  may  be  pres¬ 
ent.  Although  the  errors  due  to  background  subtraction,  splicing, 
and  uncorrected  nonlinearities  in  the  detector  responses  may  be  as 
large  as  20%,  they  are  often  smaller  and  would  be  too  subtle  to  re¬ 
cognize  without  the  normalization.  Note  that  the  corrections  for 
these  effects  described  herein  are  applied  to  the  un-normalized  data. 

The  large-scale  deviations  in  the  overall  shape  of  an  SWS 
spectmm  are  corrected  using  the  assumed  continuum  function 
as  a  guide.  Then  the  spectmm  and/or  the  12  individual  SWS 
spectral  segments  are  brought  into  agreement  with  the  accurate 
infrared  photometry  by  iteratively  removing  artifacts  and  mak¬ 
ing  finer  scale  adjustments  in  the  continuum.  Figures  1  and  2 
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Fig.  3. — Adjustments  to  a  Tau  in  the  SWS  spectral  segment  2C,  which 
contains  the  SiO  fundamental.  The  discontinuity  at  12  /rm  between  SWS  bands 
2C  and  3A  is  likely  due  to  the  nonlinear  response  of  the  detectors  or  an  over¬ 
subtraction  of  the  background.  Two  potential  corrections  for  the  discontinuity 
are  derived  by  (1)  adding  a  corrective  background  constant  or  (2)  multiplying 
the  spectral  fragment  by  2°'^,  normalized  such  that  the  correction  factor  is  1 .0  at 
7.3  /rm.  The  two  corrections  are  indistinguishable  and  overlap  as  the  gray  dotted 
lines  in  the  figure.  (This  discontinuity  does  not  appear  in  Fig.  1ft,  as  that  shows 
the  band-to-band  normalized  spectmm,  whereas  for  this  work  we  start  with  the 
non-band-normalized  spectral  segments  [Sloan  et  al.  2003].) 

graphically  demonstrate  the  before  and  after  results  for  a  Tau, 
the  first  standard  for  which  Cohen  et  al.  created  a  composite 
spectrum  (Paper  II). 

Figure  la  shows  the  Cohen  et  al.  q  Tau  composite  spectrum 
after  it  has  been  divided  by  an  Engelke  function  with  an  effec¬ 
tive  temperature  of  3898  K  (the  effective  temperature  adopted  in 
Paper  II  for  this  star),  Figure  \b  shows  the  similarly  normalized 
SWS  spectrum  from  Sloan  et  al.  (2003),  and  Figure  Ic  shows  our 
final  corrected  spectrum.  The  near-infrared  portion  of  the  Cohen 
et  al.  spectrum  (2  <  4  /rm)  is  higher  than  the  Engelke  reference, 
a  common  feature  in  many  of  the  Cohen  et  al.  spectra.  Similarly, 
the  long-wavelength  portions  of  the  normalized  SWS  spectra  are 
also  often  lower  than  the  Engelke  reference  (~20%  lower  in  the 
a  Tau  example). 

Figure  2  shows  the  analogous  comparisons  between  the  mea¬ 
sured  photometry  and  that  calculated  from  the  spectra  shown  in 
Figure  1.  Significant  deviations  between  the  measured  and  pre¬ 
dicted  photometry  occur  for  the  Cohen  et  al.  composite  spectrum 
(Fig.  2a)  at  wavelengths  below  5  /rm.  The  photometric  compar¬ 
ison  using  the  Sloan  et  al.  spectrum  (Fig.  2b)  shows  it  is  in  good 
agreement  with  the  measured  photometry  in  the  near-infrared  but 
is  discordant  by  nearly  20%  at  wavelengths  greater  than  ~1 0  /rm, 
which  is  consistent  with  the  large  uncertainty  that  Leech  et  al. 
(2003)  assign  to  the  SWS  calibration  at  the  longer  wavelengths. 
Finally,  Figure  2c  shows  that  the  deviations  of  the  measurements 
from  the  predictions  in  our  final  calibrated  spectrum  have  been 
reduced  by  using  the  photometry  to  adjust  the  SWS  spectral 
segments.  The  processing  steps  used  to  create  the  final  spectrum 
are  summarized  here  and  described  in  detail  in  the  following 
subsections. 

Some  of  the  normalized  SWS  spectra,  such  as  that  of  a  Tau  in 
Figure  \b,  show  systematic  trends  with  respect  to  wavelength. 
These  trends  often  occur  near  the  transitions  between  the  12  spec¬ 
tral  fragments.  Figure  3  shows  the  adjustment  that  was  necessary 


for  one  particular  SWS  segment  of  a  Tau  (band  2C)  that  lies 
within  the  SiO  fundamental  band.  The  solid  line  shows  the  seg¬ 
ment  after  the  Sloan  et  al.  initial  normalization.  The  end  of  the 
segment  is  about  20%  below  where  continuity  arguments  with 
the  continuum  approximation  say  it  should  be.  This  defieit  is  the 
reason  why  the  flux  from  the  longer  wavelength  speetrum  of  this 
star  in  Figure  lb  is  systematically  low.  The  droop  in  this  spectral 
segment  likely  is  due  to  ineomplete  correction  of  the  nonlinear  de¬ 
tector  response,  possibly  compounded  by  oversubtraction  of  the 
background.  The  deviation  ean  be  corrected  by  either  (1)  adding  a 
constant  to  the  background  (70  Jy  in  this  case)  or  (2)  multiplying 
by  a  weak  function  of  wavelength,  here  (2/7.3)®'^.  The  values  of 
the  offset  or  the  exponent  in  the  2"  expression  vary  with  the  wave¬ 
length  range  and  the  star  under  consideration,  but  the  eorrections 
that  we  derive  from  the  two  approaches  are  essentially  the  same. 
The  exponent  n  is  generally  small;  the  largest  value  used  for  any 
spectral  segment  for  any  of  the  stars  is  ~0.4.  We  chose  to  use  the 
2”  form  of  the  correction  for  the  general  method  of  reshaping  in¬ 
dividual  segments  for  its  ease  of  application  and  because  it  ean  ac¬ 
count  naturally  for  problems  produced  by  variations  in  response. 
The  dotted  lines  in  Figure  3  show  the  results  of  applying  the  cor¬ 
rections  to  the  a  Tau  spectral  segment;  the  two  (superposed)  cor¬ 
rections  are  indistinguishable. 

Residual  trends  and  offsets  with  respect  to  the  normalized  con¬ 
tinuum  are  corrected  in  a  piecewise  fashion  for  a  given  spectral 
segment  or  group  of  consecutive  segments.  The  corrected  stellar 
spectrum  was  thereby  initially  made  to  conform,  to  within  a  few 
percent,  to  the  Engelke  function  near  4,  7,  and  1 1  ^m,  that  is,  in 
regions  outside  of  the  well-defined  absorption  bands. 

Because  the  spectral  energy  distributions  for  stars  fall  steeply 
with  wavelength,  the  SWS  spectra  become  quite  noisy  at  the  long¬ 
est  wavelengths,  even  for  the  brightest  objects.  The  spectra  beyond 
~10  fim  are  smoothed  and  resampled  at  a  gradually  increasing 
wavelength  interval  in  order  to  preserve  a  respectable  S/N,  albeit 
at  the  expense  of  a  somewhat  degraded  spectral  resolution.  We 
chose  to  use  a  wavelength-dependent  grid  spacing  for  the  final 
spectra,  with  a  wavelength  interval  given  by  (0.0001/1  /2m)2^, 
which  provides  for  10,000  samples  over  the  1-35  (im  spectral 
range.  The  step  size  at  2.36  fim  equals  0.02%  of  the  wavelength 
(2/A2  ~  4400),  and  the  steps  are  0.34%  of  the  wavelength  at 
35  fim  (2/A2  ~  300).  Depending  on  the  observing  speed,  the  in¬ 
herent  SWS  resolution  is  2/A2  ~  500—2000  at  2.4  iim  and 
2/A2  ~  200-500  at  35  /xm  (Leech  et  al.  2003). 

The  spectral  range  over  which  SWS  data  are  used  is  indicated 
by  the  designation  S[2]  in  column  (5)  of  Tables  1  and  2.  The  SWS 
coverage  begins  at  2.36  /xm  and  ends  at  the  wavelength  specified 
within  the  brackets.  The  notation  “CVF”  indicates  that  data  from 
the  ISOCAM  CVF  are  used  for  5  Dra  and  (3  Gem.  Where  smooth¬ 
ing  does  not  sufficiently  reduce  the  noise,  a  predicted  spectral  en¬ 
ergy  distribution  derived  with  the  autoshape  templating  procedure 
(described  in  §  4)  is  substituted,  indicated  by  “A[35]”  in  that  col¬ 
umn.  Cohen  et  al.  analogously  substituted  an  Engelke  function  for 
the  noisy  or  missing  long-wavelength  portion  of  their  spectral  com¬ 
posites  and  templates. 

3.2.  Splicing  the  Near-Infrared  Data  to  the  SWS  Spectra 

The  long-wavelength  side  of  the  K  band  and  DIRBE  band  2 
just  overlap  with  the  2.36  /xm  lower  limit  to  the  SWS  spectra, 
while  the  DERBE  band  1  analog  to  the  7 band  and  the  Hammersley 
et  al.  (1998)  7 and  //photometric  data  lie  at  wavelengths  much 
less  than  the  SWS  limit.  In  order  to  make  use  of  the  near-infrared 
photometry  to  adjust  the  absolute  flux  levels,  the  calibrated  spec¬ 
tra  are  extended  to  wavelengths  shorter  than  the  SWS  limit  by 
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Fig.  4. — Calibrated  SWS  spectra  (solid  lines)  for  12  stars  ordered  by  decreasing  effective  temperature  and  normalized  by  eq,  (1)  (i.e.,  the  Engelke  function)  with 
the  appropriate  temperature  and  angular  size.  The  dashed  lines  show  the  continuum  estimates  from  eq.  (2).  The  dotted  lines  show  the  autoshape  spectmm  for  the 
star.  The  photometric  uncertainties  (1  a)  are  denoted  by  the  vertical  error  bars  for  each  measured  value,  and  the  horizontal  bars  indicate  the  >50%  response  of  the 
filter. 


seaming  the  low  spectral  resolution  (/l/AA  ~  70)  1.22-2.4  pm  ab¬ 
solute  spectra  that  Strecker  et  al.  (1979)  obtained  from  the  NASA 
Kuiper  Airborne  Observatory  (KAO)  and  Lear  Jet  Observatory. 
The  eight  stars  in  Tables  1  and  2  with  such  spectra  are  indicated 
by  “SEW”  in  column  (5)  (Strecker  etal.  also  obtained  spectra  for 
Vega  and  Sirius,  which  are  not  included,  as  we  adopt  model  spec¬ 
tra  for  these  stars).  Prior  to  splicing,  the  Strecker  et  al.  spectra  are 
first  renormalized  to  the  adopted  flux  of  Sirius  from  Price  et  al. 
(2004).  Higher  resolution  (A/AA  ~  3000)  2.02-2.41  pm  spectra 
from  Wallace  &  Hinkle  (1997)  for  q  Tau,  7  Dra,  and  p  UMa  are 
impressed  on  the  Strecker  et  al.  spectra  for  these  three  stars, 
while  the  same  is  done  for  a  Boo  with  the  data  from  Hinkle  et  al. 
(1995).  For  those  stars  without  Strecker  et  al.  data,  a  pseudo- 
spectral  energy  distribution  based  on  the  Strecker  et  al.  data  and 
near-infrared  photometry  are  used  to  extrapolate  below  the  SWS 
spectrum  to  shorter  wavelengths. 

3.3.  Photometric  Scaling 

The  rationalized  spectral  energy  distribution  is  then  iteratively 
adjusted  to  bring  it  into  agreement  with  high-quality  photometry 
on  the  star.  We  have  two  resources  for  accurate  photometry  that 
were  not  available  to  Cohen  et  al.  when  they  created  the  calibra¬ 
tion  scale  and  network.  The  first  is  the  MSX photometry  (Price 
et  al.  2004),  which  provides  a  much  needed  longer  wavelength 
lever  arm  for  the  absolute  flux  calibration  of  the  spectra.  The  sec¬ 
ond  is  the  photometry  in  the  DIRBE  Point  Source  Catalog  (Smith 
et  al.  2004).  Cohen  (1998,  Paper  IX)  confirmed  that  the  DIRBE 
photometry  is  entirely  consistent  with  their  established  zero- 
magnitude  flux  and  used  DIRBE  photometry  to  add  tertiary  stan¬ 


dards  to  their  network.  In  decreasing  order,  the  preference  for 
photometry  with  which  we  scale  the  spectra  is  MSX,  DIRBE, 
Hammersley  et  al.  (1998),  Selby  et  al.  (1988),  and,  for  the  sec¬ 
ondary  standards  in  Table  2,  the  photometry  used  by  Cohen  et  al. 
(e.g.,  that  in  Appendix  A  of  Walker  &  Cohen  1 998).  Although  the 
DIRBE  data  set  is  both  exo-atmospheric  and  extensive,  its  useful¬ 
ness  for  the  present  analysis  varies  widely  across  the  sample  of 
stars.  Since  the  DIRBE  detectors  were  large,  42'  instantaneous 
fields  of  view,  the  stellar  measurements  often  contain  background 
emission  that  degrades  the  quality  of  data  for  the  fainter  stars.  In 
addition,  contributions  from  nearby  stars  can  only  be  separated  out 
for  measurements  in  specific  scan  geometries.  Smith  et  al.  (2004) 
rejected  contaminated  scans  from  the  photometry,  but  this  fre¬ 
quently  produced  too  few  scans  at  the  longer  wavelengths  to  pro¬ 
vide  good  photometry.  For  these  reasons,  we  only  use  the  DIRBE 
near-infiared  photometry  from  bands  1-4  (i.e.,  2  <  5  gm). 

The  initial  absolute  spectral  energy  distribution  of  a  star  is 
weighted  by  the  spectral  response  of  the  filters  used  for  the  pho¬ 
tometric  observation,  and  the  results  are  integrated  over  wave¬ 
length  to  obtain  a  predicted  absolute  flux  in  each  band.  These 
integrated  in-band  fluxes  are  divided  by  the  corresponding  inte¬ 
grated  in-band  fluxes  of  the  Sirius  model  spectram.  The  resulting 
ratios  between  the  star  and  Sirius  in  that  band  are  then  compared 
to  the  ratios  actually  measured  by  MSX,  DIRBE,  Hammersley 
et  al.  (1998),  etc.  (e.g..  Fig.  2  for  a  Tau). 

If  there  is  a  significant  (>3  a)  discrepancy  between  the  calcu¬ 
lated  and  measured  ratios,  the  stellar  spectrum  is  modified  to 
bring  the  photometric  predictions  into  better  agreement  with  a 
given  measurement.  These  discrepancies  may  be  small  in  absolute 
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terms.  For  example,  the  uncertainty  in  the  relative  AT&T photom¬ 
etry  we  use  is,  for  the  most  part,  less  than  1%  and  in  many  cases 
less  than  0.5%.  The  types  of  adjustments  applied  are  as  follows: 

1 .  If  the  trends  show  a  systematic  wavelength  dependence  in¬ 
dicative  of  an  incorrect  effective  temperature  assumed  for  the 
star,  then  the  effective  temperature  in  the  normalization  (i.e.,  the 
Engelke  function)  is  adjusted. 

2.  If  the  discrepancies  show  a  uniform  bias  independent  of 
wavelength,  then  the  angular  size  is  modified. 

3.  If  the  discrepancy  is  limited  to  a  specific  range  of  wave¬ 
lengths,  then  the  spectral  segment(s)  in  that  region  is  (are)  raised 
or  lowered  as  the  photometry  indicates.  The  scale  factor  that  is 
used  to  rationalize  the  overlap  between  the  segments  being  pieced 
together  has  a  rather  large  uncertainty  because  of  the  poor  respon- 
sivity  at  the  edges  of  the  SWS  segments,  and  adjacent  segments 
are  often  rich  in  stellar  lines,  particularly  at  the  shorter  wave¬ 
lengths.  Since  the  segments  also  often  have  slightly  different  reso¬ 
lutions,  the  “best”  scaling  factor  to  combine  the  overlap  regions  is 
somewhat  ambiguous.  The  2”  correction  applied  to  a  segment  can 
also  be  changed  slightly  to  include  relatively  small  local  correc¬ 
tions  if  these  are  within  the  measurement  uncertainty. 

4.  If  the  star  is  bright  enough  to  have  both  a  well-measured 
spectrum  and  high-quality  photometry,  the  departure  of  the  spec¬ 
trum  from  the  analytic  approximation  to  the  continuum  is  ad¬ 
justed  such  that  the  low-frequency  changes  in  the  slope  of  the 
spectral  continuum  with  wavelength  [A(continuum)/A2]  are  the 
smallest  necessary  to  provide  a  good  fit  to  the  photometry.  A 
modified  continuum,  discussed  in  §  4,  is  used  for  the  fainter  stars 
that  have  less  complete  photometry. 

Since  the  measured  and  predicted  photometry  for  each,  of  the 
stars  are  expressed  relative  to  Sirius,  the  absolute  calibration  of 


the  final  spectrum  is  directly  tied  to  the  absolute  spectral  energy 
distribution  of  Sirius.  If  the  detailed  Sirius  model  that  we  adopt 
from  Paper  I  (modulo  1.01)  were  revised,  then  each  final  spec¬ 
trum  in  this  analysis  could  easily  be  revised  as  well. 

Figure  4  shows  the  absolute  spectra  for  the  12  stars  with  the 
most  accurate  photometry  and/or  the  most  extensive  SWS  wave¬ 
length  coverage.  The  spectral  types  range  from  G2  to  M5  and  are 
displayed  in  order  of  decreasing  effective  temperature  or  later 
spectral  type,  where  is  that  determined  from  our  fitting  and 
the  spectral  type  is  from  Heras  et  al.  (2002)  and  references  therein. 
Overlain  on  each  plot  is  the  photometry  with  error  bars  used  to 
adjust  the  spectra.  The  horizontal  bars  indicate  the  bandwidth  of 
the  filters  for  each  photometric  measurement. 

4.  THE  CONTINUUM  AND  SPECTRAL  TEMPLATES 

The  continua  of  the  calibrated  spectra  of  the  stars  in  Figure  4 
were  interactively  adjusted  to  follow  the  absolute  photometry. 
There  are  subtle  systematic  differences,  increasing  with  decreas¬ 
ing  temperature,  between  the  final  ealibrated  stellar  spectra  shown 
in  the  figure  and  the  corresponding  Engelke  fimction.  For  exam¬ 
ple,  the  spectra  in  later  type  stars  do  not  rise  back  to  the  reference 
level  defined  by  the  Engelke  function  (y  =  1 ;  Fig.  4,  solid  hori¬ 
zontal  line)  between  the  molecular  absorption  bands,  such  as  at 
7  fim.  On  the  other  hand,  the  continuum  increasingly  rises  above 
the  reference  level  at  shorter  wavelengths  as  the  effective  temper¬ 
ature  is  decreased. 

From  the  12  stars  in  Figure  4,  we  have  derived  an  empirical 
correction  to  the  Engelke  function  for  2  >  2.36  iim  that  is  a  bet¬ 
ter  approximation  to  the  continua.  The  shape  of  the  continuum  in 
the  modified  function  is  defined  by  a  single  input  parameter,  the 
effective  temperature  T^g.  The  significant  infrared  features  in  the 
SWS  spectra  that  arise  from  the  molecular  absorption  by  CO  and 
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TABLE  3 

Templated  Standard  Stars 


e‘  Spectra" 

Star  Spectral  Type  (K)  (mas)  [to  X  (/im)]  Photometry' 

(1)  (2)  (3)  (4)  (5)  (6) 


Q  Aur .  G4:  III:  5450  10.23  A[35]  DIR,  IRAS 

oHya .  K3II-1II  4150  9.59  A[35]  DIR.  B,  C, /A4S 

7Aql .  K3  II  4050  7.29  A[35]  DIR,  H,  S,  T 

aTrA .  K2  II  4000  9.81  A[35]  DIR,  C, /R/IS 

eCar .  K3  III  +  B2  V  3300  ±  10%  14.59  ±  10%  A[35]  DIR,  C, //MS',  2M 


“  The  uncertainty  in  T^g  is  ±7%,  and  that  in  6  is  ±4%.  Note  that  the  uncertainties  in  T^g  and  8  are  not  independent.  The  overall  absolute 
flux  uncertainty  is  usually  close  to  the  uncertainty  cited  for  8. 

"  Spectral  segment.  (A)  autoshape  (see  text). 

'  Photometry  sources.  (DIR)  Smith  et  al.  (2004);  (H)  Hammersley  et  al.  (1998);  (S)  Selby  et  al.  (1988);  (IRAS)  Beichman  et  al.  (1988); 
(B)  Bouchet  et  al.  (1989,  1991);  (C)  Carter  (1990,  1993);  (T)  Tokunaga  (1984);  (2M)  2MASS  (Stmtskie  et  al.  2006). 


SiO  are  also  found  to  be  well  parameterized  as  functions  of  the 
effective  temperature.  We  have  created  a  deterministic  procedure 
called  autoshape,  described  in  the  following  sections,  based  on 
the  correlation  of  such  key  features  with  effective  temperature 
from  the  averages  of  the  individually  calibrated  spectra. 

Autoshape  can  be  used  to  create  spectral  templates  for  a  contin¬ 
uum  of  effective  temperatures  and  has  been  effectively  validated 
for  spectral  types  between  G2  and  M5.  Since  this  approxima¬ 
tion  uses  averaged  data  from  1 2  stars,  the  autoshape  results  are 
influenced  less  by  random  measurement  error  or  by  aberrant  fea¬ 
tures  that  might  be  introduced  by  adopting  a  single  template  to 
represent  an  entire  spectral  type.  We  substitute  the  appropriate 
autoshape  data  for  noisy  or  missing  spectral  portions  for  stars 
that  have  SWS  data  (Tables  1  and  2),  use  the  function  to  make 
templates  for  stars  with  good  photometry  but  no  SWS  spectra 
(Table  3),  and  derive  estimates  of  temperature  and  angular  diam¬ 
eter  when  good  stellar  photometry  is  available  (§  5). 

4.1.  The  Shape  of  the  Continuum  above  2.36  pm 

The  deviations  between  the  calibrated  spectra  of  the  stars  in 
Figure  4  and  the  respective  Engelke  function  continua  arise  from 
a  variety  offactors,  such  as  sources  of  opacity  other  than  H“  free- 
free  absorption.  The  Engelke  function  incorporates  an  ad  hoc 
analytic  approximation  for  the  wavelength  dependence  of  the  mea¬ 
sured  disk  integrated  brightness  temperature  of  the  Sun  (Engelke 
1992).  This  was  extrapolated  to  lower  effective  temperatures  as¬ 
suming  the  infrared  opacity  is  solely  due  to  H“  free-free  absorp¬ 
tion.  In  the  derivation  of  the  function,  it  was  also  assumed  that 
the  T(t)  relationship  for  the  star  is  approximated  by  T(t)  = 
TQ{T)(Tcg/TQ).  Carbon  &  Gingerich  (1969)  validated  this  as¬ 
sumption  over  the  temperature  range  4000-6000  K. 

Atmospheric  modeling  could,  perhaps,  provide  better  estimates 
for  the  stellar  continua  than  the  Engelke  function.  The  detailed 
models  include  a  more  exact  treatment  of  opacity  and  the  variation 
of  temperature  with  optical  depth,  as  well  as  the  gravitational  ef¬ 
fects  that  Decin  et  al.  (2004)  suggest  are  necessary  to  reproduce 
the  detailed  Spitzer  Space  Telescope  Infrared  Spectrograph  (IRS) 
calibration  spectra.  Thus,  we  examined  the  infrared  spectra  gener¬ 
ated  with  detailed  stellar  atmosphere  models  and  found  that  the 
deviations  from  the  Engelke  continuum  are  similar  to  those  ob¬ 
served  in  the  normalized,  recalibrated  SWS  spectra  in  Figure  4. 
Figure  5  shows  the  normalized  difference  between  the  flux  from  a 
Kurucz  a  Boo  model  atmosphere  at  4300  and  the  correspond¬ 


ing  Engelke  function  (the  sharp  drop  at  2  <  1.6  /rm  is  due  to  the 
inapplicability  of  the  Engelke  function  below  this  wavelength,  as 
it  marks  the  transition  from  H“  free-free  being  the  dominant  opac¬ 
ity  source  to  bound-free  absorption  dominating). 

The  model  should  include  more  of  the  applicable  physics  in 
the  stellar  atmosphere  than  the  Engelke  function.  However,  the 
Engelke  function  is  actually  in  better  agreement  with  the  avail¬ 
able  spectroscopy  and  photometry  on  a  Boo,  particularly  in  the 
near-infrared,  than  either  the  Kurucz  model  atmosphere  shown 
in  Figure  5  or  the  earlier  models  for  this  star  used  by  Cohen  et  al. 
(e.g..  Paper  VII).  The  11  other  fiducial  stars  qualitatively  show 
the  same  deviations,  both  positive  and  negative,  as  seen  in  Fig¬ 
ure  5  but  at  a  smaller  amplitude.  The  dashed  lines  shown  in  Fig¬ 
ure  4  are  the  best  scaling  of  the  deviation  in  Figure  5  that  fits  the 
photometry.  These  corrections  tend  to  increase  with  decreasing 
temperature  of  the  star.  Thus,  the  photometry  for  A  >  2  pm  can 
be  brought  into  agreement  by  simply  adding  a  scaled  version  of 
the  deviations  in  Figure  5  to  equation  (1). 

Figure  6  shows  that  an  inverse  linear  fimction  of  the  effective 
temperature  provides  a  good  fit  to  the  scaling  factors  needed 
to  adjust  the  “excess”  for  the  12  fiducial  stars  in  Figure  4.  The 
dashed  line  in  the  figure  is  the  linear  least-squares  fit  to  the  data 
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Fig.  5. — Percent  deviation  of  the  Kurucz  mode!  continuum  for  a  Boo  (Tctr  = 
4300  K)  from  the  Engelke  function  for  the  same  temperature. 
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Fig.  6. — Inverse  of  the  correction  factor  (1/frac)  to  eq.  (1)  vs.  the  effective 
temperature  of  the  star.  The  fractional  correction  is  the  amount  of  deviation  of 
the  a  Boo  continuum  in  Fig.  5  needed  for  the  continuum  estimates  in  Fig.  4 
(dashed  lines  in  that  figure).  The  linear  fit  through  the  data  is  shown  by  the  dashed 
line  0.0026  ±  0.0003  Tea  —  2.99  ±  1 .08,  which  is  used  in  eq.  (2).  (If  the  G  star  at 
~5900  K,  which  is  somewhat  of  an  outlier  and  has  a  relatively  small  correction, 
were  neglected  in  the  fit,  the  line  fit  would  be  0.00277'eff  —  8.38.) 


of  0.0026rcff  -  7.99. Forwavelengthsgreaterthan2.36/im, the 
continuum  used  in  autoshape  is  the  Engelke  funetion  plus  the 
excess  in  Figure  5  divided  by  the  scale  factor 

/(A  >  2.36  fivn)  =  Eng(A,  Tetr) 

/ ,  [Kur(A,  r„Boo)/Eng(A,  TaBoo)]  -  1 1 
\  0.00267’eff  -  7.99  J’  ^ 

where  Eng(A,7’efiF)  is  given  in  equation  (1)  and  Kur(A,7’„B(,o)  is 
the  Kuruez  model  for  a  Boo  at  TcBoo  =  4300  K.  Note  that  the 
numerator  in  the  correction  factor  is  the  same  for  all  spectra. 

The  maximum  corrections  are  small,  ~3%  or  less  for  wave¬ 
lengths  greater  than  2.5  pm.  Equation  (2)  scales  the  deviations  at 
A  >  2.36  pm  in  Figure  5  by  about  0.2  for  a  KO  star,  0.5  for  K4, 
and  1 .0  for  M4. 

4.2.  Molecular  Absorption  versus  Effective 
Temperature  above  2.36  pm 

The  molecular  absorption  features  in  the  spectra  also  show 
trends  with  effective  temperature  and,  equivalently,  spectral  type. 
Absorption  profiles  for  the  CO  and  SiO  fundamental  and  overtone 
bands  have  been  extracted  from  the  calibrated  SWS  spectra  for 
the  12  stars  in  Figure  4  after  normalizing  their  spectra  with  the 
continua  defined  by  equation  (2).  The  1 2  individual  profiles  for 
each  of  the  four  molecul  ar  bands  are  averaged  to  obtain  the  mean 
absorption  profiles  shown  in  Figure  7.  These  mean  absorption 
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Fig.  7. — Averaged  molecular  absorption  profiles,  (a)  CO  fundamental  profile  averaged  over  the  spectra  in  Fig.  4.  A  correction  has  been  applied  to  remove  the 
contribution  of  the  SiO  overtone  at  4.2  fjm  (see  text).  (3)  Average  profile  for  the  CO  first  overtone,  (c)  Average  SiO  fundamental  profile,  (d)  Average  SiO  overtone 
profile. 
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Fig.  8.~ScaIing  factors  needed  to  fit  the  profiles  in  Fig.  7  to  the  stars  from  Tables  I  and  2  as  a  fiinction  of  r,frfor  (a)  the  CO  fiindamental  profile,  (b)  the  CO  first 
overtone,  (e)  the  SiO  fundamental,  and  (d)  the  SiO  first  overtone. 


profiles  are  then  scaled  to  fit  the  actual  profiles  of  each  star  with 
an  absolutely  calibrated  SWS  spectrum  in  Tables  1  and  2.  The 
scale  factors  obtained  are  plotted  against  the  derived  effective 
temperatures  of  these  stars  in  Figure  8. 

The  CO  absorption  bands  and  the  SiO  fundamental  are  ex¬ 
tracted  from  the  normalized  spectra  in  a  straightforward  fashion. 
The  SiO  overtone,  however,  has  to  be  extracted  from  the  edge  of 
the  CO  fundamental  by  taking  differences  between  the  averaged 
spectra  of  the  relatively  warm  stars  (Figs.  4a-4e)  and  the  cooler 
stars  (Figs.  4/-4/ ).  The  higher  temperature  stars  seem  to  exhibit 
little  SiO  overtone  absorption  and  can  be  scaled  relative  to  the 
spectra  of  the  cool  stars  such  that  the  CO  fundamental  is  removed 
by  the  subtraction,  thus  leaving  the  SiO  overtone. 

The  scale  factors  for  the  CO  overtone  in  Figure  8  show  an  ap¬ 
proximately  linear  decrease  with  increasing  effective  tempera¬ 
ture,  while  those  for  the  other  molecular  bands  exhibit  more 
complex  behavior.  For  these  bands,  little  variation  in  the  rather 
small  absorption  is  observed  above  a  particular  temperature  unique 
to  each  band.  The  absorption  linearly  increases  with  decreasing 
temperature  below  that  unique  temperature,  and  then  the  features 
all  saturate  at  a  temperature  of  about  3800  K,  which  produces  an 
S-shaped  curve.  The  dotted  lines  in  Figure  8  are  analj4ic  functions 
of  Teg-  that  are  fitted  to  the  data:  a  linear  function  for  the  CO  first 
overtone  and  a  hyperbolic  tangent  function  plus  a  linear  tail  for  the 
other  three  bands.  Cohen  &  Davies  (1995,  Paper  V )  examined  the 
equivalent  width  of  the  SiO  fundamental  and  found  similar  results 
(cf.  their  Fig.  10).  This  S-shaped  dependency  is  also  apparent  in 


the  equivalent  width  analysis  of  the  molecular  bands  by  Heras 
et  al.  (2002),  as  can  be  seen  in  their  Figures  5  and  6.  Aringer  et  al. 
(1997)  provide  a  theoretical  foundation  for  this  behavior  in  the 
SiO  overtone  as  a  function  of  effective  temperature  and  gravita¬ 
tional  acceleration.  Thus,  our  results  are  in  good  qualitative  agree¬ 
ment  with  the  work  of  others. 

The  average  absorption  profiles  in  Figure  7  can  be  multiplied 
by  the  scale  factors  in  Figure  8  to  fit  the  measured  profile  for  a  star 
at  a  particular  T^ff.  A  quantitative  comparison  with  other  work 
can  be  made  by  rescaling  they-axes  in  Figure  8  to  indicate  either 
the  maximum  absorption  depth  at  the  band  head  for  each  star  or 
the  equivalent  width  of  the  molecular  absorption.  For  example,  the 
maximum  absorption  depth  of  the  molecular  features  in  the  in¬ 
dividual  stars  would  be  represented  by  the  curves  in  Figure  8  if 
the  ordinates  were  multiplied  by  the  feature  minima  in  Figure  7, 
that  is,  by  ~0.25  for  the  CO  fundamental,  0.22  for  the  CO 
overtone,  0.13  for  the  SiO  fundamental,  and  0.065  for  the  SiO 
overtone. 

To  obtain  equivalent  widths,  the  plots  in  Figure  8  must  be  mul¬ 
tiplied  by  the  integrated  area  within  each  absorption  band.  For 
quantitative  comparison  with  Heras  et  al.  (2002),  the  integration 
has  to  be  limited  to  the  narrow  limits  within  each  band  that  they 
adopted.  By  comparing  corresponding  points  on  our  Figure  7  and 
their  Figures  5  and  6,  we  initially  estimated  the  integrated  areas 
to  be  0.08, 0.01 1, 0.13,  and  0.009  fim  for  the  CO  fundamental, 
the  CO  overtone,  the  SiO  fundamental,  and  the  SiO  overtone, 
respectively. 
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Integrating  the  CO  fundamental  between  the  4.30  and  4.70  iim 
limits  used  by  Heras  et  al.  produces  an  integrated  area  of  0.08  fim, 
which  is  the  value  we  had  estimated  visually.  The  scale  factor  to 
convert  the  y-axis  of  Figure  8a  to  equivalent  widths  for  the  entire 
4.0-6.7  range  spanned  by  the  CO  fundamental  band  (after  re¬ 
moving  the  SiO  overtone)  is  0.289  /im.  The  conversion  factor  to 
integrated  area  over  the  2.38-2.45  /r.m  range  used  by  Heras  et  al. 
for  the  CO  overtone  is  0.01 3  fim,  and  the  equivalent  width  of  the 
entire  CO  overtone  band  shown  in  Figure  86  is  0.096  fitn.  How¬ 
ever,  this  width  includes  the  3.02-3.40  /im  region,  which  Heras 
et  al.  note  is  affected  by  OH  absorption.  They  do  not  derive  equiv¬ 
alent  widths  for  OH  but  note  that  it  appears  in  stars  as  early  as  KO 
and  is  prominent  in  the  M  stars.  Integrating  over  the  Heras  et  al. 
7.60-9.00  fim  limits  to  the  SiO  fundamental,  we  derive  a  factor  of 
0.137  fim  to  convert  Figure  8c  to  equivalent  width;  the  equivalent 
width  of  the  entire  band  is  0.203  fim.  Integrating  the  SiO  over¬ 
tone  over  the  4.10-4.30  fim  range  adopted  by  Heras  et  al.  gives 
0.0095  fim  and  an  equivalent  width  of  0.018  fim.  Thus,  the  pres¬ 
ent  profiles  and  the  analysis  by  Heras  et  al.  are  in  excellent  quan¬ 
titative  agreement. 

The  agreement  is  not  surprising,  as  we  analyzed  about  two- 
thirds  of  the  same  stars  used  in  the  Heras  et  al.  study.  However, 
our  spectral  rationalization  does  result  in  smaller  scatter  about 
the  trends.  Since  Heras  et  al.  (2002)  found  the  same  correlations 
as  we  but  for  a  50%  larger  sample  of  stars,  the  more  precise 
trends  that  we  derive  also  apply  to  the  larger  data  set  and,  by  in¬ 
ference  from  the  random  nature  of  the  criteria  used  to  select  these 
stars  for  SWS  observation  described  by  Kraemer  et  al.  (2002),  to 
the  general  population  of  giant  stars  with  spectral  types  G2  and 
later. 

4.3.  Spectra  below  2.36  fim 

Because  of  the  discontinuity  at  1 .6  ^m  in  the  spectra  of  stars 
cooler  than  ~6000  K  (apparent  in  Fig.  5),  a  different  normali¬ 
zation  function  is  needed  to  account  for  the  trends  in  the  1 .22- 
2.5  fim  spectral  range  of  the  Strecker  et  al.  (1979)  data.  We  use 
equation  (1)  between  1.6  and  2.36  fim  and  append  a  blackbody 
with  a  variable  brightness  temperature  for  the  1 .0- 1 .6  ^m  spec¬ 
tral  region.  The  brightness  temperature  increases  from  T*  =  Tea- 
at  1 .0  fim  up  to  the  higher  brightness  temperature  of  the  Engelke 
function  at  1.6  fim',  Tt,  follows  a  wavelength  dependence  of  the 
form  (2  —  1 .0)2  between  these  two  boundary  values.  This  func¬ 
tional  form  is  based  on  the  Strecker  et  al.  stellar  continua  (after 
adjusting  them  to  the  quality  photometry),  as  well  as  theoretical 
models  (e.g.,  those  of  Kurucz). 

The  eight  Strecker  et  al.  (1 979)  spectra  from  the  KAO  for  stars 
in  Tables  1  and  2  are  separated  into  two  groups:  three  “warm” 
stars  (/?  Dra,  a  UMa,  and  /3  Gem)  and  five  “cool”  stars  (a  Boo,  7 
And,  Q  Tau,  P  And,  and  a.  Get).  The  data  are  normalized  by  the 
extrapolation  described  above  with  the  angular  sizes  and  temper¬ 
atures  determined  from  autoshape  fits  to  the  data  at  A  >  2.36  ^m. 
Based  on  the  depths  of  the  CO  and  SiO  bands,  the  equivalent  tem¬ 
perature  of  the  average  normalized  warm  stellar  spectrum  is 
4900  K  (the  simple  average  of  the  derived  individual  tempera¬ 
tures  of  the  stars  is  also  about  4900  K).  Similarly,  the  averaged 
cool  star  spectrum  has  a  temperature  equivalent  of  4040  K.  We 
model  the  (normalized)  near-infrared  segment  of  autoshape  with 
a  linear  interpolation  in  of  the  difference  between  the  two 
average  spectra  divided  by  their  temperature  difference: 

.^IR, normalized  =  hot(A)  +  (4900  -  T^s)  ' '  ^c)qq"I  ~  (3) 


Fig.  9. — Autoshape  fiinction  for  three  different  temperatures.  The  functions 
are  normalized  to  the  appropriate  Engelke  function  at  A  >  1.6  pm  and  by  a 
blackbody  with  a  brightness  temperamre  that  varies  as  a  quadratic  function  of 
wavelength,  from  the  stellar  effective  temperature  at  1.0  pm  to  the  temperature 
of  the  appropriate  Engelke  function  at  1 .6  pm.  The  strengthening  of  the  CO  and 
SiO  fundamental  and  overtone  bands  as  the  temperature  is  decreased  is  readily 
apparent,  as  is  the  An  =  1  CN  red-band  system  between  1 .4  and  1 .6  /rm  and  the 
CO  second  overtone  at  1.6  Price  (1970)  discussed  the  temperature  depen¬ 
dence  of  these  features.  The  normalization  overestimates  the  continuum  flux  for 
the  stars  for  A  &  1 .8  ^jm,  accounting  for  the  depression  in  the  normalized  spectra 
at  the  shorter  wavelengths. 

for  1.22  pm  <  k<  2.36  pm  To  test  equation  (3),  the  result¬ 
ing  spectra  were  compared  to  our  renormalized  Strecker  et  al. 
data  with  satisfactory  results. 

4.4.  Autoshape:  Application  for  Spectral  Templates 

The  corrections  and  scaling  that  created  the  calibrated  spectra 
for  the  1 2  fiducial  stars  shown  in  Figure  4  are  strongly  constrained 
over  the  entire  mid-infrared  wavelength  range  by  the  accurate 
MSX,  DIRBE,  and  Hammersley  et  al.  (1998)  photometry.  Except 
at  the  longest  wavelengths,  this  procedure  creates  spectra  for  these 
12  stars  from  observed  data  without  appeal  to  model  constracts. 
The  Engelke  function  is  used  only  for  trend  and  continuity  pro¬ 
cessing  of  the  SWS  data.  It  is  more  difficult,  however,  to  constrain 
global  shapes  for  the  fainter  stars  to  which  Cohen  et  al.  applied 
spectral  templates.  Many  of  these  stars,  including  those  in  Table  2, 
have  only  a  few  photometric  measurements.  Also,  the  uncertainty 
in  continuum  shape  will  be  larger  for  stars  that  have  significant 
reddening.  Fortunately,  Walker  &  Cohen  (2002)  found  that  a  red¬ 
dening  correction  was  necessary  for  only  two  of  the  stars  in 
Table  2:  0.03  mag  for  6  Oph  and  0.42  mag  for  tt  Aur. 

We  use  the  trends  found  in  the  12  fiducial  stars  to  develop  the 
empirical  procedure  autoshape  so  as  to  constrain  the  revisions  to 
the  SWS  data  for  the  stars  in  Table  2  that  have  only  moderate- 
quality  mid-infrared  photometry  from  IRAS  (for  which  Walker 
et  al.  [2004]  estimate  a  6%-9%  photometric  accuracy).  Auto¬ 
shape  combines  the  continuum  defined  by  equation  (2)  with  the 
average  absorption  profiles  for  the  CO  and  SiO  fundamental  and 
overtone  bands  (Fig.  7)  scaled  to  the  effective  temperature  of  a 
given  star  using  the  correlation  shown  in  Figure  8.  Equation  (3) 
allows  us  to  extend  the  spectmm  from  2.36  down  to  1.22  pm. 
The  single  input  parameter  to  the  autoshape  algorithm  is  the  ef¬ 
fective  temperature,  T^g.  Figure  9  shows  examples  of  normalized 
spectra  produced  with  autoshape  at  three  effective  temperatures. 
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Fig.  10. — Autoshape  spectrum  for  a  Hya  {dotted  line)  compared  with  the 
spectral  template  from  Paper  X  (solid  line). 

The  spectrum  derived  from  autoshape  for  the  effective  tem¬ 
perature  of  a  given  star  is  projected  as  a  scaffold  on  which  the 
SWS  segments,  if  available,  can  be  overlain.  For  those  stars  with 
little  spread  in  their  photometry,  the  depths  of  the  SWS  absorp¬ 
tions  are  the  basis  for  estimating  the  effective  temperature.  Then, 
the  angular  size  of  the  star  is  obtained  by  multiplying  equation  (2) 
at  the  derived  temperature  by  the  normalized  spectmm  and  scaling 
the  result  to  match  the  observed  fluxes.  The  autoshape  prediction 
is  also  substituted  for  noisy  or  otherwise  unusable  SWS  spectra 
beyond  10  pm  for  many  of  the  fainter  secondary  standards  in 
Table  2.  This  is  similar  to  the  Cohen  et  al.  substitution  of  Engelke 
functions  for  noisy  or  missing  spectra  beyond  13  ^m. 

We  used  the  autoshape  routine  to  create  spectral  templates  for 
a  Aur  and  7  Aql  to  test  the  efficacy  of  the  procedure.  As  a  his¬ 
torical  note.  Low  (1973)  performed  one  of  the  first  direct  mid- 
inffared  stellar  calibrations  against  a  Aur,  as  well  as  a  Tau  and 
0  And.  Autoshape  templates  were  also  created  for  three  Cohen 
et  al.  secondary  standards  that  do  not  have  SWS  measurements, 
Q  Hya,  a  TrA,  and  e  Cen,  in  order  to  provide  improved  spectral 
resolution  for  all  the  original  Cohen  et  al.  secondary  standards. 
Figure  1 0  compares  the  results  for  q  Hya  to  the  Cohen  et  al.  spec¬ 
tral  template  for  this  star.  The  autoshape  routine  can  be  used  to 
create  higher  resolution  and  more  accurate  templates  for  all  of 
the  stars  in  the  infrared  stellar  calibration  network  (Paper  X). 
Table  3  lists  stars  for  which  we  created  spectral  templates. 


As  exemplified  in  Figure  1,  Cohen  et  al.’s  spectra  and  the  tem¬ 
plates  derived  therefrom  have  systematically  higher  flux  in  the 
~l-5  pm  spectral  region  than  the  absolute  spectra  in  this  paper 
and  the  derivative  autoshape  templates.  This  discrepancy  has  also 
shown  up  in  the  recent  literature.  Cohen  et  al.  (2003,  Paper  XIII) 
created  template  spectra  for  faint  standards  to  calibrate  the  Infra¬ 
red  Array  Camera  (IRAC)  on  Spitzer.  In  their  calibration  of  IRAC, 
Reach  et  al.  (2005)  noted  a  systematic  5%-7%  difference  between 
the  camera  calibration  against  the  A  stars  compared  to  that  from 
the  K  stars;  they  decided  to  use  only  the  A  stars  for  calibration.  The 
discrepancy  disappears  when  we  repeat  the  step-by-step  analysis 
described  by  Reach  et  al.  (2005)  using  the  appropriate  autoshape 
templates.  Thus,  caution  should  be  exercised  when  using  Cohen 
et  al.’s  faint  template  spectra  for  calibration,  particularly  with  cool 
stars  (e.g.,  Ishihara  et  al.  2006). 

5.  DERIVED  PARAMETERS 

The  stellar  effective  temperature  and  angular  diameter  can  be 
derived  by  fitting  autoshape  to  the  absolutely  calibrated  spec¬ 
trum  of  the  star.  These  quantities  are  listed  in  columns  (3)  and 
(4),  respectively,  in  Tables  1-3.  The  uncertainties  are  estimated 
by  interactively  varying  the  parameters  until  the  fits,  such  as  that 
shown  in  Figure  2c,  obviously  start  to  fail.  For  many  of  the  stars 
in  Table  2,  the  overall  spectral  shape  is  poorly  constrained  by  the 
photometiy.  Therefore,  the  relationship  between  the  band  depth 
and  effective  temperature  for  the  CO  fundamental  (Fig.  8c)  is 
used  to  estimate  the  temperature.  This  typically  gives  limits  of 
±3%  for  the  temperature  range  3800  K  <  Teg-  <  5000  K  with 
somewhat  larger  uncertainties  above  and  below  these  tempera¬ 
tures.  The  corresponding  limits  on  angular  size,  which  depend  on 
the  uncertainties  in  both  the  effective  temperature  and  the  absolute 
flux,  are  then  derived. 

5.1.  Angular  Diameter 

We  compare  the  derived  angular  diameters  to  published  val¬ 
ues  obtained  by  interferometric  means  in  Table  4  and  by  indirect 
means  in  Table  5.  The  tables  list  the  number  of  stars  N  in  com¬ 
mon,  the  average  of  the  ratios  of  the  derived  to  the  published 
angular  diameters,  and  the  bias  in  the  ratio.  The  bias  is  defined  as 
(ratio  -  1  )/crrcduced.  where  cTreduccd  ”  (SmJy/N  and  o-rms  is  the  rms 
deviation  about  the  ratio  average.  Stars  of  spectral  types  G  and 
K  are  listed  separately  from  those  with  spectral  type  M,  as  a  di¬ 
chotomy  has  shown  up  between  the  two  groups  in  some  of  the 
references. 

Overall,  our  sizes  average  about  1%  larger  than  the  mean  of 
the  values  found  in  the  literature.  The  comparisons  with  pub¬ 
lished  results  on  G-  and  K-type  stars,  however,  show  much  less 


TABLE  4 

Angular  Size  Comparison:  Interferometric  Values 


Reference 

G-K  Stars 

M  Stars 

No.  of  Stars 

Ratio  ±  (Tnns 

Bias* 

No.  of  Stars 

Ratio  ±  (Tnns 

Bias’ 

Mozurkewich  et  al.  (2003) . 

10*’ 

0.999  ±  0.018 

0.2 

9 

0.979  ±  0.023 

2.7 

Di  Benedetto  &  Rabbia  ( 1 987) . 

5 

1.023  ±  0.023 

2.2 

2 

0.976  ±  0.031 

2.2 

Perrin  et  al.  (1998) . . 

2 

1.011  ±0.006 

2.6 

1 

0.981 

Dyck  et  al.  (1998;  UD) . 

6 

1.080  ±  0.070 

2.8 

5 

1.100  ±  0.022 

20 

Di  Benedetto  &  Rabbia  (1987;  UD) . 

4 

1.013  ±  0.031 

1.0 

2 

0.995  ±  0.006 

1 

Average  (excluding  UD  results) . . . 

1.012  ±  O.OII 

10 

0.979  ±  0,003 

Note, — UD  indicates  uniform  disk. 

“  Bias  =  (ratio  -  l)/(iTni,s/\/V). 

’’  Limited  to  stars  with  less  than  10%  DIRBE  variability;  see  §  6.2. 
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TABLE  5 

Angular  Size  Comparison:  Indirect  Values 


G-K  Stars  M  Stars 


Reference  No,  of  Stars  Ratio  ±  Bias*  No.  of  Stars  Ratio  ±  Bias” 


Decin  et  al.  (2000) .  7  1.006  ±  0.008  2.0  3  1.033  ±  0.026  2.2 

Blackwell  etal.  (1991) .  5  1.020  ±  0.009  5.0  3  1.030  ±  0.006  9 

Alonso  et  al,  (2000) .  6  1.010  ±  0.015  1.6  1  1.039  ±  0.026  2.6 

Manduca  etal.  (1981) .  5  1.025  ±  0.005  11  2  1.033  ±  0.027  1.7 

Di  Benedetto  (1998) .  4  0.999  ±  0.012  0.2  0 

Perrin  et  al.  (1998) .  0  ...  ...  3  1.049  ±  0.025  3.4 

Average .  1.012  ±  0.011  1.037  ±  0.008 


“  Bias  =  (ratio  -  l)/(CTnn,/A/W). 


scatter  than  the  comparisons  with  M  stars.  Possible  explanations 
for  the  increased  discrepancy  in  the  M  star  size  include  inappro¬ 
priateness  of  our  fitting  function  to  M  stars;  limitations  in  the 
methods  used  in  the  direct  measurements,  perhaps  relating  to 
limb-darkening  corrections  (discussed  below);  and  a  lack  of  con¬ 
stancy  among  M  stars  as  a  group.  A  large  disagreement  for  M  stars 
is  also  apparent  among  the  published  angular  diameters.  There  is, 
however,  some  pattern  to  the  divergence  between  our  results  and 
those  obtained  in  particular  studies  or  with  specific  methods,  such 
as  that  all  M  stars  Irom  a  specific  reference  are  either  larger  or 
smaller  than  ours.  Also,  there  are  correlations  in  the  angular  di¬ 
ameters  obtained  by  similar  methods  in  the  sense  that  M  star 
interferometric  diameters  are  larger  than  ours,  while  indirect  de¬ 
terminations  produce  smaller  M  star  diameters. 

5.1.1.  Interferometric  Sizes 

The  values  that  we  derive  for  G  and  K  stars  are  in  good  agree¬ 
ment  with  the  interferometric  diameters  from  Mozurkewich  et  al. 
(2003),  Perrin  et  al.  (1998),  and  Di  Benedetto  &  Rabbia  (1987) 
but  average  about  2%  smaller  for  the  M  stars,  particularly  for 
stars  of  spectral  type  M2  and  later.  The  Dyck  et  al.  (1998)  diam¬ 
eters  are  ~  1 0%  smaller  than  our  results,  but  their  values  are  given 
for  a  uniform  disk;  the  other  references  in  Table  4  apply  correc¬ 
tions  for  limb  darkening  using  model  expectations.  Di  Benedetto 
&  Rabbia  (1987)  published  both  uniform  disk  and  limb-darkened 
diameters.  When  their  uniform  disk  values  are  compared  to  our 
derived  values,  the  discrepancy  between  the  G  and  K  stars  and 
the  M  stars  is  greatly  reduced.  This  suggests  that  the  M  star  differ¬ 
ences  in  the  interferometer  measurements  may  be  an  artifact  of  the 
particular  limb-darkening  models  being  applied. 

There  is  empirical  support  for  the  application  of  limb-darkening 
models  to  K  and  G  stars  but  little  specific  as  to  their  use  with 
M  stars.  For  example,  Di  Benedetto  &  Bonneau  (1990)  found  that 
the  visibility  function  for  (3  And  (MO  III)  at  1 .65  and  2.2  /xm  was 
inconsistent  with  a  conventional  limb-darkened  disk.  They  pro¬ 
posed  that  the  discrepancies  could  be  explained  if  the  star  had  a 
small  bright  spot.  Haniflf  et  al.  ( 1 995)  observed  Mira  variables  and 
found  that  Gaussian  profiles  provided  a  better  approximation  to 
the  radial  brightness  profiles  than  did  uniform  disks,  although, 
besides  being  much  later  types  than  our  sample,  these  stars  also 
show  complicated  and  asymmetric  structure. 

Comparing  results  for  individual  stars:  Quirrenbach  et  al. 
(1996)  used  the  Mark  III  optical  interferometer  to  measure  the 
angular  diameter  of  q  Boo  ( K1 .5  III )  at  five  visual  wavelengths. 
They  found  that  the  measured  uniform  disk  diameters  quanti¬ 
tatively  agreed  with  wavelength-dependent  predictions  derived 
fi'om  the  limb-darkening  coefficients  of  Manduca  ( 1 979).  Our  de¬ 
rived  diameter  of  21 .1  ±  0.2  mas  is  in  good  agreement  with  their 


“true  limb-darkened  diameter”  of  21.0  ±  0.2  mas.  Hutter  et  al. 
(1989)  also  used  the  Mark  III  to  measure  24  stellar  diameters, 
albeit  before  the  calibration  and  observing  techniques  were  well 
understood  (Mozurkewich  et  al.  2003).  Our  agreement  is  good 
(within  1  £r)for/3  And  (MO  III)and7Dra(K5  III),  although  their 
a  Ari  (K2  III)  is  significantly  larger  than  other  results,  including 
ours,  which  is  consistent  with  the  Mozurkewich  et  al.  (2003) 
comment  that  the  early  Mark  III  results  may  overestimate  the 
diameters. 

Hajian  et  al.  (1998)  first  confirmed  stellar  limb  darkening  for 
the  K  giants  a  Ari  and  a  Cas  (KO  III)  with  the  Naval  Prototype 
Optical  Interferometer;  subsequently,  the  group  did  so  for  sev¬ 
eral  additional  stars  (e.g.,  Nordgren  et  al.  1 999;  Wittkowski  et  al. 
2001).  The  data  confirm  the  presence  of  limb-darkened  radial 
profiles  but  do  not  yet  allow  model-independent  determination 
of  empirical  limb-darkening  profiles.  The  authors  note  that  the 
model-dependent  limb-darkening  corrections  are  several  times 
larger  than  the  formal  uncertainties  in  the  interferometric  obser¬ 
vations.  Our  diameter  of  6.89  ±  0.14  mas  for  a  Ari  is  slightly 
higher  than  the  limb-darkened  diameter  of  6.80  ±  0.07  mas  they 
obtained  but  well  within  the  uncertainties. 

Kervella  et  al.  (2003a)  used  the  Interferometer  Commission¬ 
ing  Instmment  on  the  Very  Large  Telescope  at  K  band  to  derive  a 
limb-darkened  diameter  for  a  Cen  (GO  V)  of  8.51  ±  0.02  mas, 
the  same  value  we  derive.  Kervella  et  al.  (2003b)  give  6.01  ± 
0.02  mas  for  Sirius,  also  at  K  band.  This  is  1 .5  cr  smaller  than  the 
6.04  mas  adopted  for  Sirius  by  Cohen  et  al.  in  Paper  I  and  3  a  less 
than  the  size  of  6.07  mas  inferred  by  the  revised  flux  for  this  star 
(Price  et  al.  2004)  that  we  adopt  in  the  present  analysis. 

Ridgway  et  al.  (1982)  made  multiple  lunar  occultation  ob¬ 
servations  on  Q  Tau  (K5  III)  over  a  variety  of  wavelengths  to 
constrain  limb  darkening  and  found  a  corrected  angular  size  of 
20.88  ±0.10  mas,  in  good  agreement  with  our  diameter  of 
20.75  ±  0.21  mas.  More  recently,  Richichi  &  Roccatagliate 
(2005)  analyzed  a  combined  set  of  lunar  occultation  and  long- 
baseline  interferometric  measurements  of  a  Tau  and  found  an  av¬ 
erage  limb-darkened  diameter  of  20.58  ±  0.03  mas,  again  within 
1  cr  of  our  value. 

The  discordance  among  limb-darkening  correction  factors  is 
highlighted  in  studies  that  observed  stars  over  a  range  of  spectral 
types.  Di  Benedetto  &  Rabbia  (1987)  observed  1 1  K  and  M  giants 
at  K  band  with  the  two-telescope  interferometer  I2T  in  France. 
They  used  Manduca  (1979)  limb-darkening  models  to  derive  cor¬ 
rection  factors  ranging  from  0.975  for  P  Gem  (KO  III)  to  1 .036  for 
/3  And  ( MO  III ).  On  the  other  hand,  Perrin  et  al.  (1 998)  used  a  cor¬ 
rection  factor  of  1 .035  for  all  K1 .5  through  M8  stars  that  they  stud¬ 
ied  at  K  band.  This  value  was  also  based  on  the  same  Manduca 
(1979)  models,  supplemented  by  those  of  Scholz  &  Takeda  (1987) 
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TABLE  6 

Temperature  Comparison 


G-K  Stars 

M  Stars 

Reference 

No.  of  Stars 

Ratio  ± 

Bias’ 

No.  of  Stars 

Ratio  ± 

Bias" 

Heras  et  al.  (2002) . 

18 

1.01  ±  0.03 

1.5 

12 

1.00  ±  0.03 

0.4 

De  Jager  &  Nieuwenhuijzen  (1987) . 

16 

1.00  ±  0.03 

0.0 

12 

1.05  ±  0.02 

8.5 

Ridgway  et  al.  (1980) . 

13 

1.00  ±  0.01 

0.0 

11 

1.00  ±  0.01 

0.0 

Van  Belle  et  al.  (1999) . 

12 

1.01  ±  0.03 

1.2 

11 

1.01  ±  0.02 

1.7 

Mozurkewich  et  al.  (2003) . 

10 

1.02  ±  0.03 

2.1 

9 

1.03  ±  0.03 

3.2 

Alonso  et  al.  (2000) . 

6 

1.01  ±  0.02 

1.4 

3 

1.00  ±  0.03 

0.2 

Dyck  et  al.  (1998) . 

6 

0.96  ±  0.03 

3.2 

5 

0.96  ±  0.02 

5.5 

Manduca  et  al.  (1981) . 

5 

1.00  ±  0.01 

0.3 

2 

0.99  ±  0.03 

0.4 

Blackwell  et  al.  (1991) . 

5 

1.01  ±  0.01 

1.7 

3 

1.01  ±  0.01 

1.5 

Di  Benedetto  &  Rabbia  (1987) . 

5 

1.00  ±  0.02 

0.0 

2 

1.04  ±  0.02 

2.6 

Di  Benedetto  (1998) . 

4 

1.03  ±  0.02 

3.3 

0 

Perrin  et  al.  (1998) . 

Average . 

2 

1.02  ±  0.01 

1.01  ±  0.02 

1.6 

1 

1.02 

1.01  ±  0.02 

°  Bias  =  (ratio  - 


for  the  coolest  stars.  Mozurkewich  et  al.  (2003)  conducted  per¬ 
haps  the  most  thorough  analysis  of  limb  darkening  to  date,  using 
85  stars  with  spectral  types  from  A7  V  to  M6  III.  They  examined 
the  variation  in  measured  angular  diameters  at  four  visual  wave¬ 
lengths  and  compared  their  results  with  those  from  other  groups 
at  K  band.  They  used  correction  factors  from  more  recent  models 
by  Claret  et  al.  (1995)  and  Diaz-Cordoves  et  al.  (1995)  to  derive 
limb-darkened  diameters  from  their  uniform  diameter  measure¬ 
ments.  At  800  nm  the  corrections  range  from  about  6%  at  KO  to 
almost  10%  for  M4  stars;  the  451  nm  corrections  are  about  4% 
larger.  The  /f-band  corrections  grow  from  about  a  1%  correction 
for  K  stars  to  just  over  3%  for  the  late  M  stars.  Thus,  applying  a 
single  correction  value  to  all  stars,  even  at  K  band,  could  lead  to 
discrepancies  of  about  4%. 

Besides  the  limb-darkening  issue,  other  effects  could  contri¬ 
bute  to  the  disagreement  between  our  derived  angular  sizes  and 
those  from  other  analyses.  In  their  equivalent  width  analysis,  for 
instance,  Heras  et  al.  (2002)  specify  wavelength  limits  for  water 
vapor  absorption  but  derive  no  values  for  the  corresponding 
equivalent  width  over  this  range.  They  note  that  model  atmo¬ 
spheres  calculated  by  Decin  et  al.  (2000,  2003a,  2003b,  2003c) 
indicate  that  the  water  vapor  absorption  is  broad  and  depresses 
the  continua  of  M2-M4  stars  by  2%,  and  by  5%  for  M5  stars. 
The  fact  that  we  do  not  account  for  water  vapor  absorption  in  de¬ 
veloping  the  autoshape  routine  may  contribute  to  the  K-M  dis¬ 
crepancies  in  angular  diameters  compared  to  the  interferometric 
determinations.  A  simplistic  view  of  the  situation  is  that  the  con¬ 
tinuum  flux  from  the  lower  atmosphere  predicted  by  the  Engelke 
function  is  more  or  less  uniformly  depressed  across  the  infrared  by 
absorption  due  to  water  vapor  and  opacity  sources  other  than  H“. 
This  absorption  could  be  (partially)  taken  into  account  by  the 
algorithmic  modifications  in  equation  (2),  which  depress  the  mid- 
infrared  continuum  of  M  stars  by  2%-3%.  The  solid  angle  of  the 
star  in  equation  (2)  would  then  increase  to  compensate  for  any 
additional  absorption.  Thus,  a  2%-5%  uncompensated  absorp¬ 
tion  would  produce  a  l%-2%  low  bias  in  the  derived  angular 
diameters. 

5.1.2.  Indirect  Determinations 

The  angular  diameters  based  on  the  indirect  infrared  flux  model 
(e.g.,  Blackwell  et  al.  1991)  methods  are  all  systematically  smaller 
than  ours  (Table  5).  The  values  for  the  M  stars  are  typically  ~4% 
smaller,  and  those  for  K  stars  are  ~  1  %  smaller.  A  primary  cause  of 


this  discrepancy  is  the  different  zero-point  flux  adopted  for  the 
flux  calibrator  used  in  a  particular  study.  For  example,  Manduca 
etal.  (1981)  fitted  models  to  the  Strecker  et  al.  (1979)  1.22-5  pm 
spectra  for  five  of  the  stars  in  our  analysis.  They  rescaled  the 
Strecker  et  al.  calibration,  which  was  referenced  to  the  Vega  model 
of  Schild  et  al.  (1971),  to  match  the  Dreiling  &  Bell  (1980)  model 
for  Vega.  The  resulting  diameters  average  3%  lower  than  ours. 
This  bias  is  primarily  due  to  the  different  zero-point  flux  adopted 
for  Vega  by  Manduca  etal.  (1981),  which  falls  3%-4%  below  the 
Vega  flux  adopted  by  Cohen  et  al.  (Paper  1)  in  this  wavelength 
region.  Adjusting  the  flux  reference  to  the  absolute  flux  scale  of 
Price  et  al.  (2004)  reduces  the  average  difference  between  their 
diameters  and  ours  to  ~1 .5%.  Comparing  our  diameters  to  those 
derived  by  Cohen  et  al.,  the  ratio  for  28  stars  from  Tables  1-3  is 
1 .0 1 2  ±  0.026;  they  estimated  1  %  internal  uncertainty,  indicating 
good  agreement  despite  the  1  %  flux  rescaling. 

The  accuracy  of  the  results  derived  by  the  indirect  methods, 
including  the  current  work,  is  directly  proportional  to  the  accu¬ 
racy  of  the  adopted  absolute  flux  calibration.  The  differences  in 
the  M  star  results  obtained  from  some  indirect  methods  may 
reflect  the  calibration  uncertainty  of  their  photometry  for  these 
stars.  Alternatively,  we  may  overestimate  the  mid-infrared  flux 
for  the  M  stars  (although  neglect  of  the  water  vapor  absorption 
mentioned  above  causes  our  diameters  to  change  in  the  wrong 
direction). 

5.2.  Effective  Temperature 

We  can  also  compare  the  effective  temperatures  derived  from 
autoshape  with  those  in  the  literature  (Table  6).  The  agreement 
between  the  effective  temperatures  that  we  derive  and  those  de¬ 
termined  by  others  is  generally  quite  good.  For  a  few  studies,  the 
M  stars  show  a  larger  disagreement  than  the  earlier  spectral  types. 
Specifically,  the  comparisons  with  Mozurkewich  et  al.  (2003)  and 
de  Jager  &  Nieuwenhuijzen  (1987)  have  effective  temperatures 
that  are  3%-5%  smaller  than  ours  for  the  M  stars,  although  the 
G-K  stars  match  well. 

Ramirez  &  Melendez  (2005)  recently  used  the  infrared  flux 
method  to  derive  temperatures  and  angular  diameters  of  giant 
and  dwarf  stars.  They  too  found  that  the  direct  methods  (specif¬ 
ically  that  of  Mozurkewich  et  al.  2003)  produce  lower  tempera¬ 
tures  for  the  M  stars  compared  to  the  infrared  flux  method  results, 
and  they  suggest  that  the  limb-darkening  correction  for  the  inter¬ 
ferometry  may  be  the  source  of  this  discrepancy. 
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Fig.  1 1 . — Effective  temperatures  vs.  spectral  index  derived  in  the  present  work 
compared  to  those  found  in  the  literature.  The  spectral  index  ranges  from  -8  for 
spectral  type  G2  to  - 1  for  G9,  0  for  KO  to  5  for  K5,  and  6  for  MO  to  12  for  M6. 

De  Jager  &  Nieuwenhuijzen  (1987)  created  an  interpolation 
table  using  effective  temperature  and  spectral  types  for  268  stars 
compiled  from  the  literature.  They  included  tables  of  tempera¬ 
ture  versus  spectral  type  for  luminosity  classes  I,  III,  and  V  and 
derived  smoothly  varying  numerical  parameters  fitted  to  the 
discrete  classifications.  Table  6  compares  the  temperatures  they 
predicted  for  the  spectral  classes  present  in  Tables  1  and  2  with 
the  temperatures  derived  here.  The  average  ratio  for  K  giants  is 
1.00  ±  0.03.  As  can  be  seen  in  Figure  11,  the  temperatures  for 
the  M  giants  in  Tables  1  and  2  are  systematically  higher  with 
respect  to  the  de  Jager  &  Nieuwenhuijzen  values.  The  average 
ratio  to  M  giants  alone  is  1.05  ±  0.02,  which  corresponds  to  a 
difference  of  about  200  K.  The  discontinuity  in  the  de  Jager  & 
Nieuwenhuijzen  data  in  Figure  1 1  is  due  to  their  inclusion  of  spec¬ 
tral  types  K7  and  K9,  which  are  not  often  used  (e.g.,  Ridgway 
et  al.  1 980;  Dyck  et  al.  1 996),  and  probably  contributes  to  the  tem¬ 
perature  discrepancy  with  our  results. 

Figure  11  shows  our  derived  temperatures  (Tables  1  and  2) 
plotted  against  the  spectral  type  together  with  mean  effective 
temperatures  from  the  literature.  The  scatter  in  the  comparisons 
is  somewhat  better  for  the  spectral  range  G-K  than  for  the  M  stars, 
but  the  agreement  is  consistent  for  the  different  spectral  types, 
especially  if  the  discordant  de  Jager  &  Nieuwenhuijzen  (1987) 
values  are  excluded.  We  conclude  that  equation  (2)  is  a  good  rep¬ 
resentation  of  the  absolute  emergent  infrared  spectral  distribution 
of  late-type  stars  and  that  its  defining  parameter,  is  the  charac¬ 
teristic  temperature  of  the  star. 

6.  UNCERTAINTIES  AND  CAVEATS 
6.1.  Estimated  Uncertainties 

Price  et  al.  (2004)  derived  a  formal  rms  uncertainty  of  1 . 1  %  in 
the  direct  calibration  of  the  absolute  flux  for  a  CMa,  the  primary 
standard  star  we  use  in  this  analysis.  In  addition  to  this  uncer¬ 
tainty  are  our  uncertainty  estimates  in  seaming  the  spectral  frag¬ 
ments  together  and  in  the  photometry  used  to  scale  the  results  to 
absolute  fluxes. 

In  general,  we  estimate  the  uncertainty  in  the  local  shape  to  be 
about  5%  near  1 .22  yum  for  the  autoshape  extrapolations,  2.5% 
at  2.36  pm,  1%  between  3  and  12  pm,  and  ~2%  at  longer  wave¬ 
lengths.  These  uncertainties  are  estimated  from  comparisons  of 


Fig.  12. — Averaged  ims  scatter  of  the  deviations  of  the  final  calibrated  nor¬ 
malized  spectra  of  a  Boo,  a  Tau,  7  Dra,  /3  And,  and  7  Cm  between  2  and  35  trm 
from  the  predicted  autoshape  curves.  The  deviations  are  less  than  1  %  over  most  of 
the  wavelength  range. 

multiple  SWS  spectra  on  the  same  star,  from  comparisons  of  the 
spectra  from  different  stars  of  same  spectral  type,  and  from  the 
variations  in  the  photometry  used  in  this  analysis  that  were  ob¬ 
tained  from  different  sources.  Stars  with  several  SWS  observa¬ 
tions  typically  require  different  trend  corrections.  The  resulting 
profiles  agree  well  overall,  but  there  is  a  l%-2%  local  disagree¬ 
ment  on  the  long-wavelength  side  of  the  SiO  bands.  Discrep¬ 
ancies  between  the  normalized  spectra  for  two  stars  of  a  given 
spectral  type  usually  are  a  percent  or  less  and  are  rarely  as  large 
as  2%.  Variation  in  the  spectral  features  within  a  spectral  class 
has  also  recently  been  observed  in  Spitzer  IRS  calibration  data 
(Sloan  et  al.  2004).  For  stars  with  Strecker  et  al.  (1979)  data,  we 
estimate  the  uncertainty  to  be  ~2%  in  the  1.22-2.36  ^m  range. 

We  (partially)  account  for  water  absorption  by  the  manner  in 
which  we  tie  the  autoshape  function  to  the  photometry.  How¬ 
ever,  absorption  by  the  V2  bending  mode  at  A  =  6.55—6.70  pm 
is  apparent  on  the  wings  of  the  CO  fundamental  for  the  cooler 
stars  in  Figure  4.  There  is  some  ambiguity  as  to  how  much  water 
vapor  absorption  is  present  that  may  lead  to  the  bias  in  the  de¬ 
rived  angular  diameters  for  the  cool  stars.  A  reasonable  flux  bias 
error  is  ~2%-4%,  which  is  about  half  that  needed  to  rationalize 
the  angular  diameters. 

Figure  1 2  shows  the  rms  of  the  deviations  between  the  final  ab¬ 
solute  spectra  of  a  Boo  (K1 .5  III),  a  Tau  (K5  III),  7  Dra  (K5  Ifl), 
f}  And  (MO  III),  and  7  Cru  (M4  III)  and  the  autoshape  results. 
The  differences  between  the  autoshape  spectra  and  the  observed, 
absolutely  calibrated  spectra  are  less  than  1  %  for  most  of  the  wave¬ 
length  range. 

6.2.  Variability 

Several  of  the  stars  in  this  study  are  known  optical  variables: 
7  Dra,  Peg,  a  Cet,  /?  Gru,  and  GZ  Peg,  for  example.  Except  for 
7  Dra,  which  has  a  spectral  type  of  K5  III,  these  stars  have  spec¬ 
tral  types  of  M2  III  or  later.  As  such,  they  would  not  have  ful¬ 
filled  one  of  the  Cohen  et  al.  original  selection  criteria:  that  the 
secondary  standard  must  be  of  spectral  type  MO  III  or  earlier  to 
avoid  the  high  probability  that  the  source  is  variable  (Paper  X). 
Cohen  et  al.  relaxed  this  criterion  to  M4  III  in  order  to  include 
7  Cru  because  some  systems  required  a  calibration  source  that 
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TABLE  7 

DIRBE  Variability 


DIRBE  Band 


12  3  4 


Star 

Spectral  Type 

(%) 

(%) 

(%) 

(%) 

P  Gem . 

.  KO  III 

1 

± 

3 

1 

± 

1 

3 

± 

3 

1 

± 

2 

Q  Boo . 

.  K1.5  III 

4 

± 

3 

1 

± 

1 

4 

± 

2 

1 

± 

1 

p  UMi . 

.  K4  III 

4 

± 

3 

1 

± 

1 

3 

± 

3 

2 

± 

2 

Q  Tau . 

.  K5  III 

0 

± 

1 

2 

± 

2 

0 

± 

1 

7  Dra . 

.  K5  III 

5 

± 

3 

2 

± 

1 

3 

± 

2 

2 

± 

2 

Al  UMa . 

.  MO  III 

3 

± 

3 

1 

± 

1 

3 

± 

3 

3 

± 

3 

p  And . 

.  MO  III 

2 

± 

3 

1 

± 

1 

a  Get . 

.  M2  III 

3 

± 

3 

1 

± 

1 

3 

± 

3 

3 

± 

2 

7  Cm . 

.  M4  III 

4 

± 

3 

0 

± 

1 

3 

± 

2 

2 

± 

1 

q'  Cen . 

G2  V 

5 

± 

3 

1 

± 

6 

P  Dra . 

.  G2  II 

5 

± 

4 

4 

± 

3 

9 

± 

7 

10 

± 

10 

6  Dra . 

.  G9  III 

2 

dt 

4 

10 

± 

3 

41 

± 

7 

17 

± 

6 

S  Eri . 

.  KO  IV 

2 

± 

5 

3 

± 

6 

11 

± 

9 

16 

± 

16 

ff  Cen . 

.  KO  III 

1 

± 

2 

7 

± 

5 

6 

± 

5 

a  UMa . 

.  KO  Ilia 

5 

± 

3 

1 

± 

1 

5 

± 

2 

3 

± 

3 

^  Dra . 

.  K2  in 

4 

± 

5 

10 

± 

2 

11 

± 

4 

24 

± 

8 

a  Ari . 

.  K2  III 

3 

± 

3 

1 

± 

1 

5 

± 

3 

4 

± 

3 

7  And . 

.  K3  Ilb 

2 

± 

3 

2 

± 

1 

3 

± 

3 

3 

± 

3 

«  Tuc . 

.  K3  III 

3 

± 

3 

2 

± 

2 

3 

± 

4 

4 

± 

4 

2  Gm . 

.  K3  III 

6 

± 

5 

4 

5 

7 

± 

8 

13 

± 

17 

(T  Oph . 

.  K3  II 

4 

± 

4 

5 

± 

4 

10 

± 

7 

10 

± 

12 

i  Psc . 

.  K5  III 

8 

± 

4 

6 

± 

3 

7 

5 

14 

± 

10 

7  Phe . 

.  K4/5  III 

3 

i 

3 

6 

± 

4 

H  Sco . 

.  K5  III 

19 

± 

11 

10 

± 

13 

«  Oph . 

.  Ml  III 

1 

± 

1 

4 

± 

3 

I 

± 

2 

AE  Cet . 

.  Ml  III 

3 

± 

3 

3 

± 

5 

6  Vir . 

.  M3  III 

3 

± 

3 

1 

± 

1 

2 

± 

2 

2 

± 

2 

p  Per . 

.  M4  II 

1 

± 

1 

5 

± 

2 

4 

± 

2 

TT  Aur . 

.  M3  II 

12 

± 

6 

5 

± 

3 

P  Peg . 

.  M2.5  III 

5 

± 

3 

3 

± 

1 

4 

± 

2 

4 

± 

1 

P  Gm . 

.  M5  III 

6 

± 

3 

2 

± 

1 

4 

± 

2 

2 

± 

1 

GZ  Peg . 

.  M4  III  +  A2  V 

2 

± 

2 

6 

± 

5 

6^  Lyr . 

.  M4  III 

2 

± 

2 

2 

± 

1 

4 

± 

3 

3 

± 

2 

a  Aur ; . 

.  G4;  III: 

1 

± 

3 

0 

± 

1 

2 

± 

2 

2 

± 

1 

Q  Hya . 

.  K3  II-III 

2 

± 

3 

1 

± 

1 

2 

± 

3 

2 

± 

2 

a  TrA . 

.  K2  II 

0 

± 

2 

6 

± 

3 

5 

± 

2 

£  Car . 

.  K3  III 

3 

± 

3 

1 

± 

1 

3 

± 

3 

4 

± 

2 

7  Aql . 

.  K3  11 

4 

± 

4 

0 

± 

2 

3 

± 

4 

5 

± 

4 

Note. — DIRBE  bands  are  ( 1)1 .25  pm,  (2)  2.2  pm,  (3)  3.5  pm,  and  (4)  4.9  pm. 

was  brighter  than  any  of  those  in  the  original  network.  The  var- 

iation  in  the  infrared  brightness  of  a  star  that  is  an  optical  variable 
or  suspected  variable  is  usually  much  less  than  that  in  the  visual. 
Thus,  only  sources  with  small-amplitude  (suspected)  variability 
in  the  visual  are  included  in  Tables  1  and  2.  The  estimated  infra¬ 
red  variation  should  be  included  in  the  total  uncertainty  assigned 
to  the  calibrated  flux. 

MSX  measured  a  systematic  8%  variation  in  the  infrared  irra- 
diance  of  P  Peg  during  the  1 0  month  mission  in  1 996- 1 997.  The 
amplitude  and  phase  of  the  variation  were  well  correlated  across 
the  wavelength  range  of  4.3-21  /xm.  Smith  et  al.  (2004)  also 
statistically  detected  variability  in  this  star  during  the  10  month 
DIRBE  mission  (which  did  not  overlap  with  the  10  month  MSX 
mission).  The  0  Peg  moderate-resolution  spectrum  created  un¬ 
der  the  present  work  is  anchored  to  the  median  values  of  the  MSX 
fluxes.  To  be  conservative,  a  ±4%  uncertainty  should  be  root-sum- 
squared  with  the  ~  1  %  uncertainty  in  generating  the  spectmm  and 
the  1.1%  uncertainty  in  the  absolute  flux  calibration  and  the  result 
used  for  this  star. 


None  of  the  other  stars  in  Table  1  with  MSX  photometry  were 
measured  by  MSX  to  be  variable.  As  part  of  their  processing  to 
create  the  DIRBE  point-source  catalog,  Smith  et  al.  (2004)  sta¬ 
tistically  assessed  the  variability  of  the  bright  infrared  stars  over 
the  10  month  COBE  mission.  Table  7  lists  the  percent  variability 
in  the  DIRBE  photometry  for  bands  1  -4  from  this  catalog  for  the 
stars  in  Tables  1-3,  along  with  the  uncertainty  in  the  estimate. 

Smith  et  al.  tried  to  eliminate  scans  contaminated  by  extraneous 
sources  in  the  beam  from  .the  analysis  used  to  derived  the  variabil¬ 
ity  parameters.  However,  as  noted  previously  for  the  DIRBE  mid- 
infrared  photometry,  the  statistics  for  variability  are  less  secure  for 
the  fainter  stars.  Given  these  difiSculties  we  consider  a  source  to 
be  variable  if  the  DIRBE  variability  parameter  has  a  significance 
of  2  ff  or  greater  in  more  than  one  DIRBE  band.  Two  stars  in  the 
list,  S  Dra  and  f  Dra,  have  a  large  variability,  greater  than  10%. 
Neither  star  has  MSX  photometry.  However,  as  mentioned  earlier, 
S  Dra  was  used  to  calibrate  the  ISO  CVF.  If  the  DIRBE  results  are 
correct,  then  the  CVF  calibration  uncertainties  rnay  be  underes¬ 
timated.  Besides  0  Peg,  four  other  stars  fulfill  the  2  a  variability  in 
more  than  one  band:  a  TrA,  tt  Aur,  p  Per,  and  0  Gm.  The  DIRBE 
variability  of  0  Gm  is  somewhat  smaller  than  for  0  Peg  and  is  well 
detected  due  to  the  high-S/N  measurements  that  DIRBE  obtained 
for  this  star. 

6.3.  Outliers 

The  derived  temperatures  for  e  Car  and  6^  Lyr  are  anoma¬ 
lously  low  for  their  spectral  types;  e  Car  (K3  III)  is  a  double  star, 
but  any  contamination  from  the  hotter  B2  V  should  bias  the  de¬ 
rived  temperature  high.  Interstellar  reddening  would  bias  the  tem¬ 
perature  low,  as  would  very  thin  circumstellar  emission,  although 
Cohen  et  al.  use  zero  reddening  in  Paper  X  for  e  Car.  Regardless  of 
the  cause  of  this  anomaly,  the  final  calibrated  spectmm  of  e  Car  fits 
the  infrared  photometry  to  <  5  %.  The  source  6^  Lyr  ( M4 II )  is  also 
cooler  than  the  other  M4  stars  in  our  sample  (Table  2).  It  is  rela¬ 
tively  faint,  however,  and  has  correspondingly  higher  uncertainty 
in  the  temperature  estimate;  thus,  it  is  only  1  a  below  the  expected 
value  of  Teff  ~  3500-3600  K. 

7.  SUMMARY 

We  have  improved  the  spectral  resolution  and  absolute  accu¬ 
racy  for  10  of  the  Cohen  et  al.  secondary  standards  and  replaced 
the  q’  Cen  model  and  the  templated  spectmm  of  0  UMa  with 
absolutely  calibrated  observed  spectra.  The  moderate-resolution 
ISO  SWS  spectral  fragments  for  these  stars  from  Sloan  et  al. 
(2003)  were  scaled  to  accurate  infrared  radiometry  from  MSX 
and  DIRBE.  Low-resolution  near-infrared  spectra  from  Strecker 
et  al.  (1979)  for  nine  stars  were  seamed  onto  the  SWS  spectra  af¬ 
ter  correcting  the  data  set  to  the  Price  et  al.  (2004)  absolute  flux 
for  Sirius,  the  primary  standard  used  in  this  study.  We  adopted 
the  absolute  fluxes  and  model  spectra  proposed  by  Cohen  et  al.  in 
Paper  I  for  Vega  and  Sirius  (scaled  by  1 .0 1  for  Sirius)  to  preserve 
the  zero-magnitude  fluxes  of  Cohen  et  al.  as  validated  by  Price 
et  al.  (2004)  and  because  the  SWS  spectra  for  these  stars  become 
noisy  at  2  S;  7  ^m. 

We  have  created  a  routine  to  match  the  continuum  for  the  12 
standard  stars  as  defined  by  the  infrared  photometry.  This  routine, 
autoshape,  is  a  combination  of  the  Engelke  function  of  the  ap¬ 
propriate  effective  temperature  and  the  excess  of  a  Kumcz  a  Boo 
atmospheric  model  spectmm.  The  amount  of  excess  with  respect 
to  the  Engelke  function  is  well  fitted  by  an  inverse  linear  function 
in  effective  temperature.  The  absorption  profiles  for  the  CO  and 
SiO  fundamental  and  overtone  bands  are  found  to  vary  smoothly 
as  a  function  of  effective  temperature.  An  analytic  function  was 
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Fig.  13.2.  Final  photometric  fit  for  a  Boo. 

Fig.  Set  13. — Final  photometric  fits  for  the  stars  in  the  sample.  [See  the 
electronic  edition  of  the  Astronomical  Journal  for  Figs.  13.1-13.38.] 


derived  that,  when  combined  with  the  autoshape  continuum, 
produces  a  representative  infrared  spectrum  for  any  given  spec¬ 
tral  type.  A  pseudocontinuum  was  created  for  wavelengths  shorter 
than  the  2.36  pm  lower  limit  of  the  SWS  based  on  the  Strecker 
et  al.  (1979)  spectra  to  extend  the  autoshape  predictions  to  1 .2  pm 
for  the  stars  without  Strecker  et  al.  data.  The  autoshape  routine 
was  used  to  generate  spectral  templates  for  three  of  the  original 
Cohen  et  al.  secondary  standards  that  do  not  have  SWS  spectra  (a 
Hya,  a  TrA,  and  e  Car).  A  template  for  the  fourth  such  secondary 
standard,  p  Gem,  was  created  by  combining  the  SWS  and  CVF 
spectra  of  6  Cen,  6  Eri,  and  6  Dra,  which  are  of  the  same  spectral 
class  as  P  Gem,  with  the  autoshape  template  and  scaling  the  re¬ 
sults  to  the  high-quality  MSX,  DIRBE,  and  Hammersley  et  al. 
(1998)  photometry. 

Absolute  infrared  spectra  for  an  additional  2 1  fainter  second¬ 
ary  standards  were  created  by  combining  their  SWS  spectra  with 
autoshape  predictions  over  the  spectral  region  in  which  the  SW S 
data  were  too  noisy  to  yield  usefbl  information  and  with  the  best 
available  photometry.  Of  these  fainter  secondary  standards  with 
SWS  data,  only  a  Ari  also  had  MSX  photometry.  All  of  the  stars 
have  high-quality  DIRBE  measurements  in  at  least  two  DIRBE 
bands,  although  seven  stars  were  noted  by  Smith  et  al.  (2004)  to 
vary  by  more  than  1 0%  in  at  least  one  of  the  bands.  The  spectra 
generated  here  for  all  33  stars  will  be  made  available  via  CDS. 

The  estimated  1  a  absolute  uncertainties  in  the  secondary  stan¬ 
dards  vary  with  wavelength.  They  are  generally  l%-2%  but  can 
be  as  high  as  about  5%  at  ~1  /rm  for  spectra  that  rely  on  auto¬ 
shape.  The  minimum  uncertainty  is  dominated  by  the  absolute 
accuracy  in  the  MSW photometry  (Price  et  al.  2004).  Variability  in 
some  of  the  stars  in  Tables  1  and  2  has  been  noted,  and  the  ampli¬ 
tudes  of  the  variability  have  to  be  folded  into  the  flux  uncertainty 
estimates.  The  uncertainties  in  the  spectra  of  the  fainter  secondary 
standards  are  higher  and  reflect  the  paucity  of  high-quality  infra¬ 
red  photometry  available  for  these  stars. 

The  effective  temperatures  and  angular  diameters  of  the  stars 
were  derived  by  the  infrared  flux  method  using  autoshape  for  the 
stellar  source  fonction.  The  good  agreement  between  the  effec¬ 
tive  temperatures  that  we  derive  and  those  in  the  literature  attests 
to  the  fact  that  autoshape  provides  good  approximations  to  the 


relative  infrared  spectral  energy  distributions,  even  for  cool  M 
giants.  However,  the  opacity  sources  that  become  significant  for 
the  cool  stars,  such  as  water  vapor  absorption,  may  bias  the  de¬ 
rived  angular  diameters  to  being  too  low. 


We  wish  to  thank  the  referee  R.  Gehrz  for  helpful  comments 
that  improved  the  manuscript.  This  research  has  made  use  of 
NASA’s  Astrophysics  Data  System  Bibliographic  Services,  data 
products  from  the  Two  Micron  All  Sky  Survey  (2MASS),  which 
is  a  joint  project  of  the  University  of  Massachusetts  and  the  In¬ 
frared  Processing  and  Analysis  Center/Caltech  funded  by  NASA 
and  the  NSF,  and  the  SIMBAD  database,  operated  at  CDS, 
Strasbourg,  France. 

APPENDIX 

PHOTOMETRIC  SCALING 
OF  THE  CALIBRATED  SPECTRA 

This  appendix  contains  information  on  the  adjustments  that 
were  used  to  rationalize  the  spectra  and  to  scale  the  result  to  ab¬ 
solute  values  for  each  of  the  stars  in  Tables  1-3.  The  stars  are 
listed  in  the  order  that  they  appear  in  the  tables.  The  adjustments, 
photometric  scaling,  and  derived  parameters  are  listed  for  each 
star.  Plots  for  each  star  (Fig.  Set  13)  display  the  measured  pho¬ 
tometry  against  what  is  predicted  from  the  absolute  spectmm.  As 
an  example,  the  parameters,  photometry  (Table  8),  and  plot  for  a 
Boo  (Fig.  13.2)  are  included  here;  the  rest  are  available  in  the 
electronic  version  of  the  Astronomical  Journal. 

The  ISO  SWS  spectra  we  used  were  the  averaged  spectral 
fragments  given  by  Sloan  et  al.  (2003),  and  the  listed  adjustments 
are  applied  to  those  fragments  as  described  in  the  main  text.  The 
short-wavelength  spectra  from  Strecker  et  al.  (1979)  are  seamed 
onto  the  SWS  with  the  listed  adjustments.  Additional  short  spec¬ 
tral  fragments  from  other  resources  were  used  for  a  small  number 
of  stars  as  noted. 

For  convenience,  the  photometric  magnitude  scale  for  a  given 
photometric  reference  was  adjusted  so  that  [IRsirius]  =  -L36. 
The  DIRBE  photometry  was  taken  from  Smith  et  al.  (2004),  but 


TABLE  8 

a  Boo  Photometry 


Band 

Magnitude 

Reference 

DIRBE  1 . 

-2.221  ±  0.030  (0.04  ±  0.03) 

Smith  et  al.  (2004) 

DIRBE  2 . 

-2.988  ±  0.009  (O.OI  ±  0.01) 

Smith  et  al.  (2004) 

DIRBE  3 . 

-3.083  ±  0.022  (0.04  ±  0.02) 

Smith  et  al.  (2004) 

DIRBE  4 . 

-2.934  ±  0.020  (0.01  ±  0.01) 

Smith  et  al.  (2004) 

J . 

-2.198  ±  0.007 

Hammersley  et  al.  (1998) 

H. . 

-2.853  ±  0.007 

Hammersley  et  al.  (1998) 

K . 

-2.964  ±  0.014 

Hammersley  et  al.  (1998) 

L . 

-3.097  ±  0.010 

Hammersley  et  al.  (1998) 

J. . 

-2.22  ±  0.02 

Selby  etal.  (1988) 

K . 

-3.05  ±  0.02 

Selby  et  al.  (1988) 

L . 

-3.14  ±  0.02 

Selby  et  al.  (1988) 

MSX  A . 

-3.140  ±  0.004 

Price  et  al.  (2004) 

MSX  B, . 

-3.082  ±  0.007 

Price  et  al.  (2004) 

MSX  B2 . 

-3.03 1  ±  0.006 

Price  et  al.  (2004) 

MSXC . 

-3.187  ±  0.003 

Price  et  al.  (2004) 

MSXO . 

-3.182  ±  0.002 

Price  et  al.  (2004) 

MSX^ . 

-3.199  ±  0.008 

Price  et  al.  (2004) 

Notes. — Photometry  scaled  assuming  Sirius  is  -1.36(0.07).  Uncertainties  are 
I  <T  from  the  reference  and  taking  the  Sirius  scaling  into  account.  Numbers  in  pa¬ 
rentheses  for  DIRBE  are  the  variability  amplitudes  and  associated  uncertainties. 
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0.045  mag  was  added  to  DIRBE  band  1  results  to  correct  a  bias  for 
Sirius  in  the  DIRBE  data.  If  a  photometric  reference  had  no  Sirius 
measurement,  Vega  was  then  used,  with  magnitudes  being  calcu¬ 
lated  using  the  appropriate  bandpass  and  the  Kurucz  model  from 
Paper  I.  The  MSX  photometry  is  from  Price  et  al.  (2004),  while 
the  IRAS  information  was  obtained  from  the  VizieR  Web  site.® 
The  stellar  effective  temperature  (Teg)  and  angular  diameter  {9) 
in  milliarcseconds  were  derived  by  fitting  the  autoshape  function 
to  the  photometry.  The  plots  show  how  well  the  autoshape  func¬ 
tion  using  the  derived  parameters  fits  the  observations. 

Star: 

a  Boo 

HD  124897,  HR  5340,  IRAS  14133+1925 
K1.5  III 

CWW  composite  star 

Fit  Parameters:  T^e  =  4350  K,  6  =  21.06  mas 
Spectral  Data: 

SWS  target  dedicated  time  (TDT):  45200101  2.4-21  ^m 
Hinkle  et  al.  (1995):  2. 1-2.4  ^m 
Strecker  et  al.  (1979):  1.22-2.4  //m 


Uncertainty: 

Spectral  shape:  5%  near  1  ^m,  3%  near  2  pm,  1%  for  2.5- 
12  pm,  1.5%  beyond 
Normalization  relative  to  Sirius:  0.15% 

Spectral  Adjustments: 

SWS:  lA  0.94(2/2.5)®  ’® 

IB  0.95(A/2.6)®02 
ID  0.96(2/3.0)®  '^ 

IE  0.975(2/3.7)®  ’’ 

2A  1.06(2/4.05)®’® 

2B  1.12(2/5.35)®’® 

2C  0.975(2/7.4/  ’®  (7.0  pm  <  2  <  8.5  pm) 

2C  1.00(2/8.5)®  '®  (8.5  pm  <  2  <  10.5  pm) 

2C  1.03(2/10.4)-®-®  (10.5  pm  <  2  <  11.9  pm) 

2C0.920X  (2  >11.9  pm) 

3A  1.074(2/12.0?’ 

3C  1.06x 
3D  1.04(2/18)®-’® 

Strecker  et  al.  (1979)  0.98(2/2.0)®-®® 

Hinkle  et  al.  (1995)  1.005(2/2.2)"®-®® 


*  See  http://vizier.hia.nrc.ca/viz-bin/VizieR. 
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